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(4)  INTRODUCTION 


The  subject  of  our  project  is  the  diagnostic  imaging  of  breast  tumors  with  new  imaging 
technology,  Laser  Optoacoustic  Imaging  System  (LOIS).  The  main  merit  of  LOIS  is  significantly 
enhanced  tissue  contrast  based  on  optical  absorption.  We  are  testing  hypothesis  that  prevailing 
absorption  of  short  laser  pulses  by  advanced  angiogenesis  in  tumors  can  convert  these  tumors 
into  bright  acoustic  sources  on  the  background  of  normal  tissues  with  low  optical  absorption  in  the 
near-infrared  spectral  range.  Ultrasonic  waves,  generated  in  tumors,  can  deliver  diagnostic 
information  to  the  breast  surface,  where  signals  may  be  detected  by  piezoelectric  transducers. 
Reconstruction  of  two-dimensional  images  of  breast  sectors  can  be  made  on  the  basis  of  time- 
resolved  detection  and  analysis  of  ultrasound  signals  detected  with  an  array  of  transducers  located 
at  the  breast  surface.  Optoacoustic  tomography  takes  advantage  of  modern  wide-band 
ultrasonography  and  solid-state  laser  technology.  Particularly,  we  develop  signal  processing  and 
image  reconstruction  algorithms  based  on  methods  employed  in  wide-band  ultrasonography. 
Further  advancement  of  ultrasonic  detectors  into  ultrawide-band  detection  range  (20kHz  to  20MHz) 
opens  possibility  for  further  enhancement  of  the  optoacoustic  image  resolution. 

The  major  anticipated  advantages  of  the  laser  optoacoustic  imaging  are  as  follows:  (1)  high 
optical  contrast  between  tumors  and  normal  tissues  which  results  possibly  from  enhanced  density 
of  the  microcirculation  network  within  and  around  tumors  and  enhanced  thermoelastic  expansion 
coefficient  in  cancerous  tumors;  (2)  high  sensitivity  of  detection  owed  to  efficient  opto-thermal 
generation  of  pressure  in  tumors,  and  sensitive  detection  of  resulting  ultrasonic  waves  with  novel 
piezoelectric  transducers;  (3)  high  (<1  mm)  spatial  resolution  which  results  from  pulsed  laser 
excitation  and  time-resolved  detection  of  ultrasonic  signal  profiles  with  array  detectors;  (4) 
substantial  depth  of  monitoring  (up  to  8-cm)  which  results  from  deep  penetration  of  the  near-infra- 
red  photons  in  tissues,  efficient  conversion  of  laser  energy  into  ultrasonic  waves  propagating 
through  the  breast  with  insignificant  attenuation  and  minimal  distortion,  and  from  high  sensitivity 
(>2  V/bar)  of  piezoelectric  transducers:  and  (5)  applicability  to  radiographically  dense  breasts. 


(5)  BODY  OF  REPORT 

The  purpose  of  this  BCRP  project  is  to  incorporate  laser  excitation  into  a  wide-band 
ultrasonography  system  (with  a  linear  array  probe)  to  yield  a  clinical  prototype  Laser  Optoacoustic 
Imaging  System  (LOIS)  and  test  this  system  initially  in  phantoms  and  then  in  breast  cancer 
patients.  The  tissue  phantoms  with  embedded  blood  vessels  were  designed  to  resemble  natural 
optical  and  acoustic  properties  of  breast  and  tumors  with  advanced  angiogenesis.  Patients  with 
racial  differences  in  skin  types  were  included  in  the  studies  to  account  for  variation  in  optical 
absorption  properties.  Clinical  studies  performed  during  the  Year-2  of  the  project  confirmed  our 
hypothesis  that  the  breast  structural  heterogeneity  in  vivo  does  not  compromise  the  optoacoustic 
image  quality  and  may  be  accounted  for  in  the  optoacoustic  image  processing. 

During  the  second  year  period  of  the  project  all  specific  aims  and  milestones  were  accomplished. 
LOIS  performance  was  improved  by  incorporating  multichannel  electronics  system  allowing  real¬ 
time  image  acquisition.  Initial  clinical  tests  of  the  system  were  made  on  breast  cancer  patients. 
Specific  details  of  the  optoacoustic  imaging  system  development  and  tests  are  presented  below. 

Task  3  (Year  2):  Evaluation  of  the  system  using  50  patients  with  advanced  cancer  and  scheduled 
for  radical  mastectomy  (complete  surgical  removal  of  the  breast). 

Milestones  of  Task  3  were  accomplished.  The  system  performance  was  tested  in  vivo  on  10 
patients  chosen  from  the  examination  group  of  50  patients  diagnosed  with  breast  cancer  and 
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scheduled  for  the  radical  surgical  mastectomy  Not  all  oatients  were  eligible  for  the  optoacoustic 
imaging  studies  due  to  the  fact  that  their  tumors  were  very  large,  occupied  substantial  portion  of  the 
breast,  were  palpable  and  were  often  treated  with  radio-  and  chemo-  therapy.  These  large  tumors 
were  clearly  visualized  with  x-ray  radiography  and  had  significant  radiographic  contrast.  Therefore, 
such  tumors  being  an  “easy  target”  for  conventional  breast  imaging  technologies,  could  not  be 
used  to  test  hypotheses  outlined  in  our  proposal.  The  following  hypotheses  were  tested  (i)  whether 
LOIS  is  capable  of  detecting  malignant  breast  tumors  m  vivo,  (ii)  whether  contrast  of  optoacoustic 
images  is  greater  than  the  contrast  of  x-ray  radiography  and  ultrasound,  and  (iii)  whether 
optoacoustic  images  of  malignant  tumors  may  be  differentiated  from  that  of  benign  tumors. 

Fig.  1  presents  typical  optoacoustic  images  obtained  in  vivo  from  breast  cancer  patients. 
Procedure  was  performed  in  the  surgery  suite  before  radical  mastectomy  procedure. 


Figure  1.  Optoacoustic  image  of  a  malignant  tumor  obtained  /n  vivo  (left,  UTMB  patient  #  X15221X), 
Optoacoustic  image  of  a  malignant  tumor  obtained  in  vivo  (right,  UTMB  patient  # 
X26372X).  Note  that  the  contrast  is  so  significant  that  the  entire  background  of  normal 
tissue  appears  dark. 


B 

Figure  2.  X-ray  mammography  image  of  a  cancerous  breast  obtained  in  vivo  with  General  Eelectric 
mamography  system,  Mo  anode,  Mo  filter  (A)  There  were  2  suspicious  areas  on  mammogram, 
only  one  of  them  was  a  ductal  carcinoma  Ultrasound  image  of  the  same  breast  showing 
malignant  tumor  obtained  with  Acuson  128  XP/iO  3  5  MHz  USI  system  (B).  Note  that  contrast  of  the 
tumor  relative  to  background  is  not  sufficient  to  define  exact  shape  and  dimensions  of  this 
malignancy. 


The  contrast  and  precision  of  tumor  localizaiio''!  n  optoacoustic  images  were  compared  with 
that  of  conventional  x-ray  mammography  and  ultrasound  images.  The  same  breast  was  imaged  in 
vivo  with  General  Electric  mammography  system  with  Mo  anode  and  Mo  filter  and  with  Acuson  128 
XP/10  system  employing  3.5  MHz  C3-curved  linear  array  transducer.  These  images  are  depicted  in 
Fig.  2.  X-ray  image  revealed  two  suspicious  areas  with  no  microcalcifications.  A  core-needle  biopsy 
performed  after  x-ray  screening  proved  that  only  one  of  the  two  suspicious  areas  was  cancerous, 
and  the  other  area  of  the  breast  similarly  bright  on  mammography  image  was  found  normal.  The 
ultrasound  image  of  the  same  tumor  was  not  very  impressive.  The  core  of  the  tumor  and  its 
boundaries  can  not  be  clearly  seen  in  this  image 

Optoacoustic  images  were  also  correlated  u  the  gr  jss  and  histology  /  pathology  of  the 
mastectomy  specimens.  After  optoacoustic  imaging  procedure,  a  surgical  excision  of  the  breast 
with  tumor  was  performed  (Dr.  Declan  Fleming)  and  given  to  a  pathologist  (Dr.  Zoran  Gatalica)  who 
examined  the  mastectomy  specimens  in  order  to  determine  tumor  geometry,  location  of  tumors, 
the  type  of  the  tumor  and  (in  selected  cases)  microvessel  density.  Laser  safety  of  the  system  was 
also  evaluated  by  pathologist.  The  optoacoustic  imaging  procedure  was  performed  using  laser 
fluence  of  30-50  mJ/cm^  which  corresponds  to  one-third  to  one-half  of  maximum  of  permissible 
laser  exposure  for  the  skin  in  the  near  infrared  spectral  range,  as  described  in  American  National 
Standards  for  Safe  Use  of  Lasers.  Pathology  examination  compared  10  specimen  of  tissue 
involved  in  optoacoustic  imaging  and  10  specimen  of  tissue  not  involved  in  optoacoustic  imaging. 
No  thermal  or  mechanical  damage  to  normal  tissue^s  suf’'ounding  the  tumors  were  found  in 
association  with  optoacoustic  imaging. 

In  order  to  perform  optoacoustic  imaging  in  reai  '  rne  niy  a  single  optical  fiber  was  placed 
on  the  breast  surface  in  locations  thought  to  be  above  or  near  the  tumor.  Usually,  4  optoacoustic 
images  were  acquired  from  each  breast  examined  m  vivo.  Thirty  two  optoacoustic  signals  were 
employed  for  data  collection  and  subsequent  image  reconstruction.  Optoacoustic  contrast  between 
the  tumor  and  surrounding  normal  tissues  in  vivo  exceeded  200%  in  all  cases  examined  and  was 
found  equal  to  that  earlier  measured  in  vitro  mastectomy  specimen.  This  contrast  is  so  significant 
that  surrounding  normal  tissue  look  dark  on  the  image  Such  a  contrast  permits  detection  of  small 
2-mm  tumors  at  the  depth  of  6-cm. 


Figure  3.  Optoacoustic  signal  amplitude  in  aqueous  solution  4  nanoparticles  as  a  function  of 
incident  optical  energy  (A).  Note  signal  enhancement  after  the  threshold  optical  fluence  of  10 
mJ/cm  is  reached  on  the  particles.  Optoacoustic  image  of  two  small  objects  imbedded  in  the 
tissue-like  gel  phantom  (B).  Object  #1  on  the  right  was  colored  with  absorbing  nanoparticles,  while 
object  #2  on  the  left  was  colored  with  dye  solution  of  equal  absorption  coefficient.  In  spite  of  equal 
absorption,  the  brightness  of  the  object  #1  is  10  times  higher  than  the  brightness  of  the  object  #2. 


In  the  course  of  the  project  we  invented  possibility  '  urthe?  enhance  optoacoustic  contrast  with 
gold  nanoparticles  that  may  be  targeted  selectively  to  malignant  tumors.  Thermodynamics  teaches 
that  heat  capacity  of  the  particles  is  proportional  to  the?  particle  diameter  cubed,  while  the  absorbed 
portion  of  incident  optical  energy  decreases  only  as  the  riartine  diameter  squared.  Thus,  decrease 
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in  particle  diameter  yields  gradually  increasing  temperature  of  the  particle  subjected  to  constant 
fluence  of  laser  energy.  This,  in  turn,  permits  particle  superheating  without  heating  of  surrounding 
tissue.  In  order  to  overcome  the  threshold  for  evaporation,  photodissociation  or  other  gas 
production  at  a  safe  level  of  laser  fluence,  the  absorbing  particles  should  have  dimensions  on  the 
nanometer  scale  (30-150  nm).  The  optimal  size  of  nanoparticles  from  the  standpoint  of 
optoacoustic  efficiency  is  simultaneously  just  right  for  extended  circulation  of  the  nanoparticles  in 
the  blood  stream.  Intravenous  injection  of  colloidal  gold  has  been  employed  for  treatment  of 
rheumatoidal  arthritis  for  more  than  three  decades  and  proved  to  produce  tolerable  toxicity. 
Therefore,  we  anticipate  minimal  or  no  accumulation  of  the  proposed  nanoparticulate  contrast 
agent  in  normal  cells.  Gold  is  known  a  noble  metal  causing  no  visible  toxicity  in  micromolar 
concentrations.  The  progress  of  the  proposed  research  project  will  depend  on  success  of 
development  of  biomolecule-nanoparticle  coupling  technique,  allowing  reliable  conjugation  of  our 
nanoparticles  with  ligand  peptides  and  antibodies  specific  to  target-cells.  A  new  grant  application 
was  submitted  to  BCRP-2001  requesting  financial  support  for  expansion  of  the  optoacoustic 
imaging  project  with  nanotechnology. 

In  order  to  demonstrate  the  utility  of  laser  optoacoustic  imaging  for  the  quantitative 
measurement  of  the  hemoglobin  oxygenation  in  vascularized  tumors,  we  have  performed  pilot 
experiments  with  gel  phantoms  containing  embedded  “arteries  and  veins”. 

The  main  chromophores  in  tissue  that  absorb  radiation  in  the  near-infrared  range  are  oxygenated 
or  deoxygenated  molecules  of  hemoglobin  in  the  blood.  The  contrast  in  laser  optoacoustic 
imaging  is  a  result  largely  of  the  fact  that  blood  is  unevenly  distributed  in  tissue.  Growing 
cancerous  tumors  develop  an  extensive  network  of  blood  vessels  through  the  process  of 
angiogenesis.  These  vessels  are  “leaky,”  and  blood  concentrates  in  the  tissue  of  tumor  tissue. 
When  laser  optoacoustic  imaging  is  performed  with  more  than  one  wavelength  of  light,  it  has  the 
power  to  reveal  not  only  where  blood  is  concentrated,  but  the  level  of  blood  oxygenated,  measured 
by  the  relative  concentration  of  oxygenated  and  deoxygenated  hemoglobin  molecules. 

The  pooled  blood  in  cancerous  tumors  generally  tends  to  be  oxygen  poor.  The  tumor  consumes 
abnormally  large  amounts  of  oxygen  in  the  process  of  active  growth.  Furthermore,  leakage  of  the 
blood  from  the  arteries  of  the  tumor  prevents  its  recirculation  to  the  lungs  for  oxygen  replenishment. 
Thus,  relatively  few  of  the  hemoglobin  molecules  in  the  blood  in  a  cancerous  tumor  have  an 
attached  oxygen  molecule. 

Near-infrared  radiation  with  a  wavelength  of  1064  nm  is  absorbed  about  10  times  more  strongly  by 
oxyhemoglobin  molecules  (those  with  an  attached  oxygen  molecule)  than  by  deoxyhemoglobin 
molecules  (those  free  of  oxygen).  Laser  optoacoustic  imaging  with  this  wavelength  largely 
provides  a  map  of  the  distribution  of  oxyhemoglobin. 

Near-infrared  radiation  with  a  wavelength  of  757  nm  is  absorbed  more  strongly  by 
deoxyhemoglobin  than  by  oxyhemoglobin  molecules  although  both  absorb  sufficiently  strongly  to 
contribute  to  a  laser  optoacoustic  image.  Measurements  made  at  both  1064  nm  and  757  nm 
provide  the  necessary  information  for  mapping  both  blood  concentration  and  blood  oxygenation. 

The  figure  below  illustrates  the  sensitivity  of  laser  optoacoustic  imaging  to  blood  oxygenation.  The 
phantom  consisted  of  collagen  gel  mixed  with  whole  milk  that  had  been  warmed  up  and  diluted 
with  enough  water  to  give  it  optical  properties  similar  to  what  the  human  breast  would  have,  and 
then  cooled  down  to  solidify.  Enough  gelatin  was  dissolved  to  create  a  solid  mass  with  a  diameter 
of  12  cm  comparable  to  most  human  breasts.  Imbedded  in  the  phantom  were  loops  of 
polyethylene  tubing  to  simulate  blood  vessels.  In  these  loops  was  sheep  blood  that  had  been 
exposed  to  oxygen  in  varying  degrees.  As  shown  on  the  figure,  the  loop  on  the  left  contained  blood 
that  independent  measures  showed  was  37%  oxygenated.  In  other  words,  37%  of  the  hemoglobin 
molecules  (OaHb)  had  an  attached  oxygen  molecule.  The  loop  on  the  right  was  partitioned,  and  the 
top  half  contained  28%  oxygenated  blood,  and  the  bottom  half  contained  92%  oxygenated  blood. 

The  known  absorption  properties  of  oxyhemoglobin  and  deoxyhemoglobin  indicate  that  a  laser 
optoacoustic  image  taken  with  a  wavelength  of  757  nm  will  show  most  intensely  those  portions  of 
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the  phantom  that  are  rich  in  deoxyhemoglobin.  These  are  the  portions  of  the  loops  that  are  low  in 
oxyhemoglobin.  The  laser  optoacoustic  image  on  the  right  of  the  figure  indeed  shows  the  top 
portion  of  the  loop  on  the  right  most  intensely,  the  loop  on  the  left  somewhat  less  intensely,  and  the 
bottom  portion  of  the  loop  on  the  right  least  intensely. 


Laser  Irradiation  f  757  nm 


Single  Fiber.  757  nm 


Contrast  Gradient 
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Figure  4.  Schematic  diagram  (a)  and  two-dimensional  optoacoustic  image  of  two  blood  vessel 
loops  in  a  gel  phantom  with  optical  properties  similar  to  those  of  breast  tissue  in  vivo.  Imaging  was 
performed  using  Alexandrite  laser  at  760-nm. 

The  test  image  shows  clearly  the  power  of  laser  optoacoustic  imaging  to  resolve  shapes 
and  sizes  of  blood-filled  tumors  in  the  human  breasts.  It  also  illustrates  the  potential  for 
measurement  of  the  blood  oxygenation  levels  of  such  tumors.  This  experimental  image  also 
contains  a  few  imperfections.  Most  of  the  observed  artifacts,  such  as  the  bright  bands  through  the 
left  loop  result  from  the  use  of  a  single  irradiation  point  and  a  rather  small  number  of  transducers 
in  the  detector  array.  These  will  become  negligible  when  a  newer  detector  array,  under 
construction,  is  implemented  in  the  course  of  the  Year-3  of  the  project.  The  bright  band  at  the 
bottom  of  the  image  near  the  detector  array  will  be  removed  by  a  new  filtration  procedure. 

Data  of  initial  clinical  experiments  were  analyzed,  and  recommendations  were  made  for 
LOIS  modification  in  order  to  improve  its  specifications  and  prepare  for  clinical  studies  in  the 
radiology  suite  of  patients  with  smaller  breast  tumors 


Figure  5.  New  LOIS  computer  interface  for  real-time  imaging. 
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In  order  to  optimize  optoacoustic  system  characteristics  form  the  standpoint  of  diagnostic 
efficacy,  safety,  resolution  and  sensitivity,  such  laser  irradiation  parameters  as  dynamic  range  and 
ultrasonic  bandwidth  of  the  acoustic  detection  were  varied.  Optimal  parameters  of  the  array 
transducer  were  found  to  be  as  follows:  (1)  ultrasound  detection  bandwidth:  20  kHz  to  4  MHz, 
dimensions  of  each  transducer  35  mm  x  1mm  with  1  mm  gaps,  64  piezoelements  total  in  120° 
linear  array  with  diameter  of  120  mm.  Sensitivity  of  such  transducer  made  of  PZT-5H  ceramic 
embedded  in  polymer-epoxy  resin  with  2-2  structure  was  exceptionally  high,  >2  V/bar  with  effective 
noise  pressure  of  only  4  pbar.  In  the  course  of  the  next  year  project  we  also  plan  to  use  two 
different  wavelengths  of  laser  irradiation,  so  that  diagnostic  capability  of  LOIS  can  be  tested. 


(6)  KEY  ACCOMPLISHMENTS 

•  Laser  Optoacoustic  Imaging  System  (LOIS)  was  modified  and  enhanced  to  operate  in  real  time 
and  at  two  wavelengths  of  Nd:YAG  or  Alexandrite  lasers,  which  should  enable  diagnostic  capability 
of  optoacoustic  imaging. 

•  LOIS  was  tested  in  gelatin  phantoms  simulating  dense  breast  with  blood  vessels  having  blood 
with  various  degree  of  oxygen  saturation.  Experiments  demonstrated  that  LOIS  is  capable  of  not 
only  visualizing  shape  and  dimensions  of  blood  vessels,  but  also  differentiating  deoxygenated 
blood  from  oxygenated  blood. 

•  LOIS  was  tested  in  breast  cancer  patients  scheduled  for  radical  mastectomy.  Tests  demonstrated 
that  LOIS  can  visualize  malignant  tumors  in  vivo  with  contrast  significantly  exceeding  the  contrast  of 
x-ray  mammography  or  ultrasound  images. 

(7)  REPORTABLE  OUTCOMES 

•  Manuscripts:  3  published,  2  in  press,  2  submitted,  1  in  preparation.  Abstracts  -  7,  Presentations 
acknowledged  USAMRMC  support  -  7. 

•  Patents:  1  patent  application  USPTO  #  60/147,577  entitled:  “Optoacoustic  Monitoring  of  Blood 
Oxygenation”  (with  co-inventors:  Rinat  Esenaliev,  Donald  Prough,  and  Massoud  Motamedi). 
Inventor’s  response  to  the  US  Patent  Office  action  was  submitted. 

•  Based  on  work  supported  by  this  award  the  following  grant  applications  were  submitted: 

1)  US  Civilian  Research  and  Development  Foundation,  Co-PI,  “Development  of  algorithms  and 
software  codes  for  3D  optoacoustic  tomography”,  $56,800  (funded). 

2)  Cancer  Research  Foundation  of  America,  “Quantitative  tomography  of  early  detection  of 
Colorectal  Cancer”,  $40,000  (awarded). 

3)  USAMRMC  Prostate  Cancer  Program-2001,  “Novel  imaging  system  for  early  detection  of 
prostate  cancer”,  $558,600  (pending). 

4)  USAMRMC  Breast  Cancer  Program-2001,  “Diagnostics  of  microscopic  breast  cancer  using 
smart  nanoparticles  as  an  optoacoustic  contrast  agent”,  $402,489  (pending). 

5)  NIH/NCI,  “Diagnostic  imaging  of  blood  in  breast  tumor  angiogenesis”,  $1,049,936 
(pending). 

(8)  CONLUSIONS 

•  Laser  Optoacoustic  Imaging  System  is  a  useful  modality  for  detection  of  breast  cancer,  which 
provides  enhanced  tissue  contrast.  Limits  of  small  tumor  detection  and  diagnostic  capability  of  LOIS 
may  be  established  in  the  course  the  Year-3  project. 
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ABSTRACT 

A  clinical  prototype  of  the  laser  optoacoustic  imaging  system  (LOIS)  was  employed  for  breast  cancer  detection  and 
localization  in  patients  with  confirmed  breast  cancer  and  scheduled  for  radical  mastectomy.  The  prototype  LOIS  used  a 
single  optical  fiber  for  delivery  of  laser  pulses,  an  arc  shaped  32-eiement  PVDF  transducer  array  for  ultrawide-band 
piezoelectric  detection  of  optoacoustic  signals  and  a  single-channel  data  acquisition  card  for  signal  processing.  The 
resonance  ultrasound  frequency  of  the  110  pm  PVDF  film  was  outside  detectable  range  of  ultrasound.  wSpatial  resolution  of 
the  transducer  array  was  slightly  better  than  Imm  in  radial  direction  and  slightly  worse  than  1  mm  in  lateral  direction.  The 
system  was  optimized  for  contrast  and  sensitivity.  Data  acquisition,  signal  conditioning  and  image  processing  were 
significantly  improved  and  optimized  resulting  in  reduced  image  frame  rate  of  2  seconds  employing  700  MHz  Aphlon 
processor.  The  computer  code  for  digital  signal  processing  employed  band-pass  hyper-Gaussian  filtering  and  denoising.  An 
automatic  recognition  of  the  optoacoustic  signal  detected  from  the  irradiated  surface  was  implemented  in  order  to  visualize 
the  breast  surface  and  improve  the  accuracy  of  tumor  localization.  Radial  back-projection  algorithm  was  employed  adopting 
combination  of  integration  along  spherical  wavefronts  and  integration  along  planar  wavefronts  (as  in  Radon  transform)  for 
image  reconstruction.  The  system  performance  was  evaluated  initially  in  breast  tissue-like  phantoms  with  embedded  blood 
vessels.  Clinical  studies  in  breast  cancer  patients  scheduled  for  surgical  mastectomy  were  performed  and  compared  with  x- 
ray  radiography,  ultrasound  and  pathology  reports. 

Keywords:  Optoacoustic  imaging,  breast  cancer,  ultrawide-band  acoustic  transducer,  laser,  ultrasound 


1.  INTRODUCTION 

Laser  optoacoustic  imaging  system  (LOIS)  was  proposed  for  cancer  detection  in  order  to  combine  advantages  of  optical 
contrast  and  sensitive,  high-resolution  ultrasonic  detection  in  one  tomography  system  [1-4].  One  of  the  major  applications  is 
the  detection  of  early  breast  cancer,  especially  in  younger  women.  The  optoacoustic  tomography  utilizes  analysis  of  profiles 
of  acoustic  signals  induced  by  laser  pulses  in  breast  tissue  with  tumors  [5-8].  Application  of  the  ultrawide-band  ultrasonic 
detection  instead  of  detection  of  photons  in  optoacoustic  tomography  helps  to  overcome  two  problems  associated  with  strong 
light  scattering  in  biological  tissues  and  improve  depth  of  monitoring,  sensitivity  and  spatial  resolution  [9-12].  The  profiles 
of  acoustic  waves  generated  under  irradiation  conditions  of  temporal  pressure  confinement  in  the  volume  of  tumors  resemble 
the  profile  of  absorbed  laser  energy  in  the  tumors,  and  can  be  used  to  replicate  the  tumor  structure.  Tumors  with  dimensions 
of  L5-mm  to  15-mm  irradiated  with  laser  pulses  represent  themselves  as  sources  of  acoustic  waves  with  ultrasonic 
frequencies  of  ~1  MHz  to  -100  kHz,  Such  ultrasonic  waves  can  propagate  in  biological  tissues  with  insignificant  attenuation 
[12-13].  However,  spherical  propagation  of  acoustic  waves  changes  initial  profile  of  absorbed  laser  energy  (especially  in  the 
ultrasonic  frequency  range  <  1  MHz),  so  that  the  detected  signals  must  be  processed  by  integrating  them  over  the  entire  time- 
course  of  detection  in  order  to  be  used  in  the  optoacoustic  image  reconstruction  [13-15]. 

The  first  few  years  of  our  research  in  the  area  of  optoacoustic  imaging  of  the  breast  yielded  the  following  results:  (1)  first 
clinical  prototype  system  for  breast  imaging  was  developed,  fabricated  and  extensively  tested  in  phantoms  [10,12],  (2) 
detection  and  localization  of  breast  tumors  in  surgical  mastectomy  specimens  were  demonstrated  [15],  (3)  parameters  of 
optoacoustic  were  compared  with  X-ray  mammography,  MRI  and  ultrasound  imaging  [15,16].  The  main  merit  of 
optoacoustic  imaging  was  found  to  be  the  significant  tissue  contrast.  The  image  contrast  is  defined  as  the  difference  in 
optoacoustic  image  brightness  in  tumors  and  in  adjacent  normal  tissue  divided  to  the  brightness  of  normal  tissue  background. 
This  contrast  was  found  in  the  range  of  1  to  5,  which  substantially  exceeds  any  other  endogenous  tissue  contrast  currently 
utilized  in  clinical  ultrasonography,  MRI  and  X-ray  mammography  (usually  much  less  than  1).  Based  on  literature  data  and 
our  gross  observations  of  tumor  cross-sections  we  hypothesize  that  the  optoacoustic  contrast  results  primarily  from  increased 
optical  absorption  in  the  dense  micro  vascularity  of  the  tumors  [12,17].  However,  optical  scattering  in  tumors,  especially 
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those  treated  with  chemo-  and  radiotherapy,  may  also  significantly  contribute  to  the  optoacoustic  contrast  [18]. 
Thermoelastic  properties  of  the  breast  tumors  may  also  add  to  the  optoacoustic  contrast,  however,  were  not  yet  quantified 
[19].  The  lateral  resolution  of  acoustic  transducer  array  was  improved  relative  to  that  in  the  system  used  previously  in  vitro. 
Comparison  of  optoacoustic  images  with  those  of  X-ray  radiography,  MRI  and  ultrasound  imaging  performed  in  vitro 
revealed  good  correlation  in  tumor  size  and  location,  and  the  spatial  resolution  equal  to  that  of  ultrasound  (better  than  1  mm) 
was  achieved  [15]. 

In  this  paper,  we  report  initial  results  of  in  vivo  clinical  studies  involving  patients  with  breast  cancer.  The  clinical 
optoacoustic  system  LOIS  was  assembled  of  a  compact  Nd:  YAG  laser,  fiberoptic  light  delivery  system,  an  arc-array  of  32 
wide-band  acoustic  transducers,  and  a  data  acquisition  board  operated  by  computer.  Design  of  the  new  system  is  compact 
and  suitable  for  performing  studies  in  vivo.  New  signal  and  image  processing  employed  permitted  reconstruction  of 
optoacoustic  images  in  1  second  after  data  collection. 


2.  MATERIALS  AND  METHODS 


2.1.  LOIS  for  clinical  studies 

Schematic  diagram  of  LOIS-02  is  shown  in  Fig.  ia.  A  Nd:YAG  laser  (Big  Sky  Lasers,  MA)  operating  at  the  wavelength  of 
1064  nm  was  used  as  a  source  of  near-infiared  pulses  of  10-ns  duration.  Repetition  rate  of  laser  pulses  was  20  Hz.  A  quartz 
optical  fiber  was  employed  for  the  laser  pulse  delivery  to  the  tissue  surface.  A  telescope  was  used  to  expand  laser  beam  from 
the  fiber  and  deliver  a  parallel  beam  with  diameter  of  1-cm  to  the  surface  of  skin.  A  glass  window  was  in  contact  with  tissue 
during  irradiation  assuring  the  exact  position  and  providing  rigid  surface  conditions  for  acoustic  pressure  coupling.  Rigid 
surface  of  the  irradiated  tissue  permitted  more  effective  generation  of  higher  ultrasound  frequencies  at  the  very  surface  for 
better  determination  of  the  surface  position  in  the  detected  signals.  This  method  of  irradiation  also  resulted  in  generation  of 
minimal  tensile  pressure,  yielding  safe  laser  irradiation  conditions  [6]. 

2  Arc  array  of  ultrwide-band  acoustic  transducers 

A  specially  designed  array  of  PVDF  acoustic  transducers  was  employed  m  LOIS.  The  array  had  32  rectangular 
piezoelectric  transducers  of  1.5  x  12  mm  dimensions  and  4  mm  distance  between  transducers  (see  Fig.  1).  Piezoelectric 
polymer  PVDF  of  1 10-^m  thickness  was  used  for  transducer  fabrication.  Low  acoustic  impedance  and  ability  to  operate  in 
the  ultrawide  ultrasonic  frequency  band  without  strong  resonance  and  reverberations  are  the  advantages  of  PVDF  for 
detection  of  opto-acoustic  profiles.  Backing  material  was  used  for  mechanical  matching  of  transducer  and  damping 
reverberations  after  detection.  The  transducers  were  located  on  arc  surface  of  120-mm  diameter.  A  limited  viewing  angle  of 
120  degrees  resulted  in  reduced  lateral  resolution  of  1.5-mm  compared  with  in-depth  resolution  of  0.8-mm.  On  the  other 
hand,  this  geometry  provided  sufficient  resolution  in  a  wide  area  around  the  focal  zone. 


Figure  1.  Schematic  diagram  of  clinical  LOIS  (a)  and  a  photograph  of  the  32-element  arc  transducer  array  (b). 

Absolute  sensitivity  of  each  transducer  in  the  array  was  measured.  Results  of  this  calibration  were  used  in  normalizing 
signals  to  the  sensitivity  of  each  transducer.  The  mean  sensitivity  value  for  transducer  array  of  5  =  10  p.V/Pa.  Electric 
capacity  of  each  transducer  was  approximately  =  40  pF,  which  resulted  in  effective  noise  level  of  20  pV  as  shown  in 
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Figure  2.  Electrical  voltage  corresponding  to  the  effective  noise  pressure,  Ap,  can  be  calculated  for  a  transducer  operating  in 
the  idler  (open  circuit)  mode  as  follows  [18]: 


Uno,se=Sl^P  =  ^-^^^  (1) 

where  g  is  the  piezoelectric  modulus,  1  is  the  thickness  of  the  piezoelement,  kXo  is  the  room  temperature  and  Cp^  is  the 
electric  capacity  of  the  piezoelement..  Simultaneously  this  value  of  electric  capacity  for  sufficient  for  effective  coupling  of 
electrical  signals  from  transducers  into  the  charge  preamplifier  and  signal  processing  circuitry. 


Figure  2.  The  effective  thermal  noise  pressure  in  PVDF  acoustic  transducers  as  a  function  of  transducer  electric  capacity. 

2.2.  Signal  Filtering  in  LOIS 


The  major  part  of  signal  processing  in  LOIS  invloves  various  filtering  procedures.  Each  signal  detected  by  acoustic 
transducers  in  the  array  consists  of  the  following  three  components:  (1)  a  sharp  peak  produced  at  the  very  subsurface  layer  of 
the  irradiated  tissue  (a  product  of  tissue  local  absorption  and  effective  subsurface  laser  fluence),  (2)  a  low-frequency  smooth 
exponential  slope  produced  by  the  attenuation  of  light  intensity  inside  tissue  and  (3)  short  N-shaped  pulses  generated 

by  smaH  absorbing  tumors.  The  duration  of  the  N-shaped  pulses  is  defined  by  the  time  of  sound  propagation  through  the 
dimension  of  the  tumor  in  the  direction  of  detection.  The  time-delay  of  the  optoacoustic  signal  arrival  depends  on  tumor 
location  relative  to  the  transducers  in  array.  Therefore  position  of  the  tumor  and  its  dimensions  can  be  determined  from  the 
optoacoustic  signals  detected  by  LOIS. 
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Figure  3.  Typical  optoacoustic  signal  detected  by  a  central  transducer  in  array  from  a  breast  with  tumor  and  the  same 
signal  after  integration  procedure  converting  the  signal  into  a  profile  of  absorbed  energy  distribution  (left  panel),  and 
the  optoacoustic  signal  zeroed  down  at  the  surface  position  (right  panel). 
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The  signal  generated  by  laser  pulses  in  the  subsurface  layer  (very  close  to  the  surface)  possesses  maximum 
amplitude  over  the  entire  time-period  of  detection.  Therefore,  the  signal  gradient  is  maximal  at  the  tissue  surface,  which 
yields  high  ultrasound  frequencies.  High  ultrasound  frequencies  that  fall  near  the  resonance  frequency,  of  piezo¬ 
elements  can  induce  significant  reverberation  in  the  transducer.  Our  transducers  were  designed  to  damp  resonance 
frequencies  and  widen  the  range  of  detectable  frequencies.  Nevertheless,  a  strong  reverberation  takes  place  after  the  arrival 
of  the  signal  from  the  surface  of  tissue.  These  reverberations  do  not  affect  the  useful  part  of  the  optoacoustic  profile. 
Besides,  it  is  convenient  to  employ  the  sharp  peak  with  high  gradient  (derivative)  for  automatic  determination  of  the  tissue 
surface.  The  surface  position  can  be  accurately  detected  at  the  intersection  of  integrated  signal  with  zero  line  (see  Fig.  3,  left 
panel).  The  controlling  computer  generated  numerical  signal  at  the  surface  position  determined  for  each  transducer  and 
stored  these  data  in  a  separate  file.  Ail  data  samples  were  equaled  to  zero  after  the  arrival  of  the  signal  generated  at  the  tissue 
surface.  The  surface-marks  generated  by  computer  were  used  for  the  surface  visualization.  Combined  image  contained 
surface  marks  and  the  image  of  breast  tissue  with  tumors,  showing  tumor  location  relative  to  the  irradiated  surface. 

Before  the  integrated  signal  could  be  used  for  image  reconstruction,  it  is  often  desirable  to  eliminate  or  reduce  its  low 
frequency  component  associated  with  homogeneous  light  attenuation  inside  the  tissue,  because  it  can  significantly  decrease 
the  image  contrast.  This  was  achieved  in  LOIS  with  high-pass  Gaussian  filtering  with  RC  cells.  The  value  of  RC^-IO  ps 
defines  the  cut-off  ultrasonic  frequency,  and  the  sharpness  of  the  filter  cut-off  is  defined  by  the  number  of  RC  cells.  The 
greater  the  number  of  elementary  RC  cells,  the  steeper  the  slope  of  the  filter  transfer  function.  The  cut-off  frequency  and  the 
slope  of  the  transfer  function  of  the  filter  could  be  varied  conveniently  by  the  radiologist-operator. 


Normalized  Frequency, 


Figure  4.  Hyper-Gaussian  filter  transfer  function  (left  panel)  and  filter  responses  for  various  types  of  RC  filters  (right  panel). 


In  case  when  large  tumors  are  present  in  the  breast,  the  optoacoustic  signals  coming  from  tumors  carry  significant  low- 
frequency  component.  Therefore,  elimination  of  exponential  trend  associated  with  effective  optical  attenuation  in  tissue  may 
simultaneously  remove  useful  information  about  tumor  internal  structure  from  the  image.  Ideal  filter  would  be  represented 
by  a  step- function  0-to-l  with  infinitely  sharp  slope.  However,  real  filter  possess  certain  resolution,  yielding  smooth  slope 
from  0-to-l  over  a  period  of  its  RC  constant.  It  is  important  to  note  that  the  standard  RC  filter  can  have  a  level  greater  than 
zero  only  after  it  is  initiated.  The  “inertia”  of  a  standard  RC  filter  produces  not  only  widening  of  signals  after  filtration,  but 
also  delays  of  the  signal  components  relative  to  their  original  positions  on  the  time  course.  These  delays  mean  phase  shifts  in 
the  frequency  domain  and  result  in  blurring  of  optoacoustic  images  reconstructed  using  signals  filtered  with  a  standard  RC- 
filters.  Figure  4  (left)  shows  graphically  the  difference  between  an  ideal  step-function  and  a  real  filter.  To  avoid,  the 
problem  of  phase-shifts  upon  filtering,  we  employed  Hyper-Gaussian  filters  with  variable  parameters  for  band-pass  filtering 
of  integrated  optoacoustic  signals: 


G  =  exp[(-^r]  (1) 

Jeff 

where  G  is  the  filter  transfer  function,  /  is  the  ultrasonic  frequency  in  optoacoustic  signals,  is  the  effective  cut-off 

frequency,  and  n  is  the  hyper-Gaussian  parameter.  The  graphical  expression  of  this  filter  is  depicted  in  Figure  4  (right  panel). 
Integrated  signals  filtered  from  high  frequency  noise  were  used  as  input  data  for  the  image  reconstruction  code. 


Proc.  SPIE  Vol.  4256 


9 


Application  of  the  filters  shown  in  Fig.  4  for  elimination  of  high-frequency  noise  and  low-frequency  exponential  trends  in 
optoacoustic  signals  yielded  profiles  without  widening  or  delays  in  signal  components.  Figure  5  depict  integrated 
optoacoustic  signal  before  and  after  filtering  low-frequency  components  with  hyper-Gauss ian  filter.  In  spite  of  generally 
satisfactory  results,  some  distortions  can  be  observed  after  filtration,  showing  artifacts  of  negative  signal  values,  minor  inside 
tissue  and  especially  noticeable  above  the  tissue  surface.  Any  artifacts  in  the  signals  were  equaled  to  zero  before  the  image 
reconstruction. 

Depth,  mm 

0  20  40  60  80  1  00  1  20 


Figure  5.  Integral  of  a  band-pass-filtered  optoacoustic  signal  starting  at  the  tissue  surface  and  the  same  integral  filtered 
through  a  low-pass  Gaussian  filter  in  order  to  eliminate  exponential  trend  associated  with  effective  optical  attenuation  in 
tissue  and  reveal  the  position  and  the  size  of  a  breast  tumor, 

23,  Optoacoustic  Image  Reconstruction 

A  radial  back-projection  algorithm,  developed  for  the  optoacoustic  imaging  and  described  in  our  earlier  publications  was 
employed  [10,  20-22].  For  the  reconstruction  of  two-dimensional  optoacoustic  images  we  used  signal  integrals  shown  in  Fig 
5  and  projected  them  back  onto  the  two-dimensional  grid  taking  into  account  the  directivity  pattern  (angle  of  acceptance)  of 
each  transducer  in  the  airay.  The  image  represents  distribution  of  a  product  of  thermo-acoustic  efficiency,  optical  absorption 
and  absorbed  laser  energy.  The  back-projected  images  acquired  with  only  limited  number  of  detectors  display  artifacts 
associated  with  and  image  reconstruction  using  incomplete  set  of  data.  The  contrast  of  these  images  could  be  improved 
through  the  image  filtration  procedure  apphed  to  the  entire  two-dimensional  image  [15].  However,  unprocessed  optoacoustic 
images  contain  quantitative  information  of  distribution  of  absorbed  optical  energy  in  tissue,  This  quantitative  information 
may  be  used  for  functional  imaging  and  estimate  of  tissue  optical  properties. 


Figure  6.  LOIS  computer  interface. 


LOIS  interface  is  presented  in  Fig.  6.  Each  new  image  is  displayed  on  the  left  panel.  The  nght  panel  allows  accumulation 
and  averaging  of  multiple  images.  Data  acquisition  time  in  LOIS-02  equals  13  seconds  for  one  fiber  position.  The  laser 
pulse  repetition  rate  is  assumed  to  be  20  Hz  and  16  pulses  for  averaged  acquisitions  are  employed.  Data  processing  and 
back-projected  image  formation  takes  0.8  seconds.  Image  filtration  and  its  optimal  processing  take  0,2  seconds.  Total  time 
equals  14  sec.  Current  design  of  the  laser  optoacoustic  imaging  system  iXOIS)  allows  operator  to  scan  the  acoustic  detector 
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array  manuaily  through  the  breast  and  to  obtain  images  of  chosen  parts  ot  the  breast. 

3.  RESULTS  AND  DISCUSSION 

3.1.  Initial  Tests  of  Clinical  LOIS 

Laser  optoacoustic  imaging  system  was  tested  in  various  phantoms.  As  shown  in  the  diagram  on  Figure  7,  the  artificial  blood 
vessels  were  located  at  different  depths  and  were  placed  either  parallel  or  perpendicular  to  the  imaging  plane.  Position  of  the 
transducer  array  was  at  the  bottom  of  the  image  (not  shown).  Surface  position  is  visible  as  a  bright  area  on  the  images.  A 
single  fiber  was  used  for  irradiation  at  only  one  location  on  the  surface  of  the  phantom.  A  photograph  of  the  phantom  is 
presented  in  Fig.  7  (right).  The  phantom  was  imaged  twice,  one  time  with  optical  fiber  positioned  at  one  point  of  its  surface 
and  the  second  time  at  the  opposite  position  (the  phantom  was  rotated  1 80  degrees). 


Figure  7.  Left  panel:  Collagen  gel  phantom  with  artificial  blood  vessels  filled  with  rabbit  blood  having  various  levels  of 
oxygenation.  Positions  of  blood  vessels  and  percent  of  oxyhemoglobin  are  shown  in  windows.  Optical  properties  of  this 
phantom  were  similar  to  the  optical  properties  of  the  breast  at  the  wavelength  of  1064  nm.  Right  panel:  Photograph  of  the 
same  phantom  with  artificial  blood  vessels.  Transducer  array  and  optical  fiber  are  also  presented  on  the  photograph. 


Contrast  Bar  Contrast  Bar 


Figure  8.  Optoacoustic  image  of  blood  vessels  in  milky  gelatin  phantom.  Contrast  of  blood  vessels  relative  to  background  is 
shown  on  Contrast  Bar.  Depth  from  the  surface  is  depicted  on  the  left  vertical  axis  of  each  image. 

Two-dimensional  images  of  blood  vessels  of  different  diameter  and  filled  with  blood  having  different  level  of  oxygenation  is 
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shown  in  Figure  8.  Whole  milk  was  diluted  by  water  to  obtain  the  effective  optical  attenuation  coefficient  1.4  cm“'  in  order  to 
mimic  the  upper  limit  of  optical  scattering  in  the  breast  tissue  at  the  wavelength  of  1064  nm.  Optical  absorption  coefficients 
of  blood  in  blood  vessels  varied  from  0.8  to  4.0  per  cm  depending  on  oxygen  content.  The  first  image  was  reconstructed 
using  optoacoustic  signal  integrals  filtered  through  a  low  frequency  filter  in  order  to  remove  the  exponential  trend  associated 
with  effective  optical  attenuation  (Fig.  8  left).  The  second  image  was  not  filtered  through  a  low  frequency  filter,  so  that 
exponential  profile  of  light  distribution  in  the  phantom  could  be  observed  (Tig.  8,  nght). 

The  image  of  blood  vessel  cross-sections  differs  from  the  conec:  circular  shape.  This  fact  has  a  simple  explanation.  As  the 
angle  of  acceptance  of  these  objects  by  the  transducer  array  was  about  120°  only  one-third  of  data  needed  for  complete  image 
reconstruction  were  collected  and  employed  in  image  reconstruction.  This  limited  lateral  resolution  in  the  far  zone  of  the 
images  presented  in  Fig.  8.  An  increase  of  the  acceptance  angle  for  the  blood  vessel  located  farther  from  the  irradiation  site 
and  closer  to  the  detector  array  resulted  in  a  better  reproduction  of  its  circular  cross-section  in  the  image,  as  clearly  depicted 
in  Fig.  8  (right  panel)  showing  an  image  of  a  blood  vessel  in  the  near  zone.  The  angle  of  acceptance  in  this  case  is  about  180° 
and  as  a  consequence  the  shape  of  image  highly  correlates  with  real  shape  of  the  circle.  Note  that  our  image  reconstruction 
algorithm  allows  accurate  reproduction  of  spherical  objects  and  cylindrical  objects  (see  a  blood  vessel  in  the  center  going  in 
horizontal  direction).  Thus,  we  may  predict  that  objects  with  random  shapes  will  be  correctly  reproduced  on  LOIS  images, 
The  low  frequency  filtration  procedure  demonstrated  in  Fig.  8  (left  panel";  yield  image  with  completely  eliminated  artifacts 
due  to  1/r  noise-enhancement  near  transducers. 


All  blood  vessels  were  clearly  resolved  and  depicted  as  separate  objects  The  relative  position  of  the  blood  vessels  as 
depicted  on  the  optoacoustic  images  accurately  resembles  their  position  as  shown  on  photograph  Fig.  7  (right  panel)  and  their 
depths  are  in  good  agreement  with  their  real  location.  This  imaging  expenment  demonstrated  sensitivity  of  LOIS  sufficient 
for  detection  of  blood  vessels  at  the  depth  of  over  75  mm  in  an  optical  strongly  scattering  and  absorbing  phantom.  The 
brightness  of  the  blood  vessel  images  was  proprtional  to  the  amount  of  absorbed  optical  energy,  which  in  turn  was 
proportional  to  the  level  of  blood  oxygenation  in  each  vessel 

3.2.  Quantitative  Optoacoustic  Imaging 


Optoacoustic  tomography  visualizes  the  spatial  distnbution  of  absorbed  optical  energy  in  tissue.  However,  distortion  of 
optoacoustic  signals  upon  propagation  through  thick  layers  of  tissue  and  signal  processing  may  in  pnnciple  modify  the 
original  profile.  We  tested  possibility  to  obtain  axial  in  depth  profiles  of  absorbed  optical  energy  from  optoacoustic  images 
from  human  muscle  tissue.  A  human  arm  was  imaged  in  the  area  of  byceps  with  LOIS  designed  for  breast  cancer  imaging. 


1  '  I  ■  I  •  [  ^  [ 
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Figure  9.  Optoacoustic  image  of  human  arm  in  the  area  of  byceps  ( left  panel )  and  in-aepth  cross-section  of  this  image 
resembling  the  profile  of  absorbed  optical  energy  (right  panel  ).  The  effective  optical  attenuation  coefficient,  il^^  =  1.28  cm' 7 


Onl>  high  frequency  filtenng  was  applied  to  the  optoacoustic  image  'fnis  expenments  demonstrated  that  quantitative 
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information  can  be  obtained  form  optoacoustic  images.  Neither  signai  correction  for  acoustic  diffraction,  nor  reconstruction 
of  absorbed  optical  energy  distribution  from  optoacoustic  signals,  nor  correction  for  sensitivity  of  individual  piezoelectric 
transducers  changed  quantitative  information.  Tne  effective  optical  attenuation  coefficient,  cm'\  determined  form 

exponentil  fit  of  the  optoacoustic  profile,  was  in  a  good  agreement  with  results  of  direct  measurements  of  optical  properties  in 
this  tissue. 


5.2.  Imaging  Breast  Cancer  in  vivo 


Eligible  patients  were  chosen  from  the  group  of  patients  ^cheaulcd  tor  radical  surgical  mastectomy  with  breast  cancer 
diagnosed  with  x-ray  mammography  or  a  combination  ot  x~ray  mammography  and  other  imaging  modalities  (ultrasound, 
MRl)  and  biopsy.  The  patients  were  imaged  in  the  surgery'  room  before  surgery.  Two-dimensional  optoacoustic  images 
were  acquired  in  several  locations  on  the  breast  with  cancer.  Irradiation  with  512  pulses  from  a  nanosecond  NdiYAG  laser 
operating  at  the  wavelength  of  1064  nm  was  performed  only  at  one  sue  or  breast  skin  surface.  Optoacoustic  transducer  array 
was  placed  on  the  opposite  side  of  the  breast  approximately  beneath  the  point  of  irradiation.  Irradiation  point  was  placed 
approximately  above  the  area  suspicious  of  being  a  tumor.  The  exact  location  of  tumors  was  not  known  prior  to  opto¬ 
acoustic  imaging  procedure,  however,  approximate  location  could  be  determined  from  the  mammogram.  Tumors  were 
sometimes  palpable  and  biopsy  incision  was  visible  on  the  spared  segment  ot  skin.  However,  some  tumors  were  not  palpable 
and  located  deep  within  the  breast.  Two  exemplary  optoacoustic  images  of  breast  carcinoma  is  presented  in  Figure  10. 
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Figure  10.  Optoacoustic  image  of  a  malignant  tumor  obtained  m  vivo  >.  left.  LTMB  patient  #  315221Q),  Optoacoustic  image 
of  a  malignant  tumor  obtained  in  vivo  (right,  UTMB  patient  #  026372P). 


Figure  11.  X-ray  mammography  image  of  a  cancerous  breast  (left).  There  were  2  suspicious  areas  on  mammogram,  only  one 
of  them  was  a  ductal  carcinoma.  Ultrasound  image  of  the  same  breast  showing  malignant  tumor.  UTMB  Patient  #315221. 
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The  main  conclusion  that  could  be  made  from  clinical  experiments  on  breast  cancer  patients  is  that  optoacoustic  tomography 
provides  enhanced  contrast  between  normal  tissues  and  cancerous  tumors.  This  contrast  varies,  however,  no  one  tumor  of  5 
tumors  examined  had  contrast  less  than  1,  i.e.  the  optoacoustic  amplitude  is  twice  higher  in  tumors  relative  to  normal 
background.  Optoacoustic  images  presented  in  Fig.  10  depict  details  of  tumor  structure  based  mainly  on  optical  properties  of 
cancerous  tissues.  Location  of  the  tumors  inside  the  breast  and  the  tumor  dimensions  also  were  accurately  determined  from 
the  optoacoustic  images  and  confirmed  through  comparison  with  x-ray  mammography,  ultrasound  and  pathological 
examination  (see  Figure  11).  Figure  11  presents  x-ray  mammography  and  ultrasound  images  of  the  breast  with  a  malignant 
tumor,  which  was  depicted  in  optoacoustic  image  in  Fig.  10  (left).  The  position  of  the  tumor  relative  to  the  skin  surface  was 
better  determined  from  ultrasound  image.  However,  detailed  structure  of  the  tumor  was  better  visualized  in  x-ray 
mammography  image.  On  the  other  hand,  x-ray  image  depicted  two  areas  suspicious  of  being  cancerous,  and  only  one  area 
contained  malignant  tumor  (as  determined  by  pathology  study).  The  tumor  core  brightly  displayed  in  the  optoacoustic  image 
was  not  visible  in  neither  x-ray  or  ultrasound  images,  but  was  confirmed  after  subsequent  pathology  examination.  The 
dimension  of  a  tumor  core  can  be  estimated  as  10x7  mm.  Tumor  was  detected  at  the  depth  of  11  mm  from  the  surface. 
Contrast  between  the  tumor  and  surrounding  normal  tissues  in  this  image  exceeds  1.2.  The  contrast  in  ultrasound  image  and 
x-ray  mammography  image  of  the  same  tumor  did  not  exceed  0.1.  Resolution  of  image  visualization  with  LOIS  (--  1-mm) 
was  comparable  with  that  of  x-ray  mammography  and  ultrasound.  Results  of  comparison  between  optoacoustic,  ultrasonic 
and  x-ray  radiography  images  for  the  patient  026372  (see  Fig.  10,  right)  and  for  the  other  3  patients  were  qualitatively 
similar. 


4.  CONCLUSION 

The  first  clinical  prototype  of  the  laser  optoacoustic  system  for  two-dimensional  imaging  of  the  breast  in  vivo  was  developed 
and  employed  for  detection  of  breast  cancer  in  5  patients.  The  system  employed  an  array  of  32  elements  of  ultrawide-band 
PVDF  transducers  (20kHz- 2MHz)  shaped  as  an  arc  of  120  degrees.  The  system  initially  was  tested  in  various  phantoms  and 
demonstrated  specifications  that  make  it  suitable  for  high-contrast  imaging  human  tissues.  Sensitivity  of  LOIS  permits 
detection  of  2-mm  blood  vessels  at  the  depth  of  7.5  cm.  In  depth  resolution  equals  0.4  mm,  lateral  resolution  is  about  1  mm 
depending  on  position  of  tumor  relative  to  the  transducer  array.  Quantitative  information  on  distribution  of  absorbed  optical 
energy  can  be  obtained  from  optoacoustic  images.  This  clinical  prototype  LOIS  is  the  basis  for  the  next  generations  of  three- 
dimensional  imaging  systems  to  be  developed  at  UTMB  in  collaboration  with  LaserSonix  Technologies,  Inc.  (Houston,  TX). 
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ABSTRACT 

Mathematical  model  of  image  reconstruction  for  two-dimensional  optoacoustic  imaging  system  is  described.  It  was  assumed 
that  (1)  receiving  transducers  are  uniformly  distributed  along  the  perimeter  of  a  60-mm  radius  ring  with  2.1 -mm  gaps 
between  transducer  centers  and  (2)  initial  data  were  known  with  0.1-mm  increments.  The  algorithm  of  radial  back  projection 
with  convolution  was  used  for  optoacoustic  image  reconstruction.  The  convolution  was  evaluated  with  modified  Shepp- 
Logan  (MSL)  and  rectangular  (RECT)  space  spectrum  windows.  Linear  interpolation  was  applied  for  calculation  of  the 
convolution  at  the  intermediate  space  points.  The  following  four  criteria  were  employed  for  estimation  of  resulting  image 
quality:  (1)  noise  level  on  entire  tomogram,  (2)  a  jump  transfer  function,  (3)  loss  contrast  function  and  (4)  the  contrast- 
dimension  relation.  Theoretical  expressions  for  these  parameters  were  derived  and  used  for  optimization  of  the  proposed 
algorithm.  Two-dimensional  images  of  computer  simulated  spherical  objects  were  reconstructed.  It  was  shown  that  0.1-mm 
spatial  resolution  could  be  obtained  provided  the  signal- to-noise  ratio  equals  approximately  3  at  the  tomogram.  A  very  small 
(0.2-mm  diameter)  tumor  and  a  small  (2 -mm  diameter)  tumor  could  be  clearly  revealed  at  the  tomogram  if  their  optical 
absorption  contrast  equals  at  least  2  and  0.1  respectively. 

Keywords:  image  reconstruction,  tomography,  filtered  backprojection  algorithm,  image  quality  criteria 

1.  INTRODUCTION 

A  clinical  prototype  of  the  laser  optoacoustic  imaging  system  (LOIS)  with  32-arc  transducer  array  was  successfully 
tested  in  mastectomy  samples  and  in  vivo  studies  [1,2].  Pronounced  optoacoustic  contrast  of  1.2-to  3.5  was 
demonstrated.  Design  of  the  acoustic  array  allows  operator  to  scan  through  the  breast  manually  and  obtain  images  of 
chosen  parts  of  the  breast  in  real  time.  The  main  drawback  of  this  system  is  a  limited  viewing  angle  (120^),  which 
resulted  in  insufficient  lateral  resolution.  A  new  optoacoustic  tomography  system  designed  for  stationery  investigation 
is  under  development  in  UTMB  now.  It  employs  acoustic  array  with  180  transducers  uniformly  distributed  along  the 
surface  of  hollow  cylinder.  A  breast  is  placed  inside  cylinder  filled  with  degassed  and  deionized  water.  It  is  illuminated 
by  laser  pulses  delivered  through  optical  fibers  located  in  gaps  between  acoustic  transducers.  A  complete  set  of  data  for 
high  quality  image  reconstruction  can  be  collected  in  this  configuration.  A  purpose  of  this  work  was  to  develop  an 
algorithm  and  software  code  for  the  breast  cancer  image  reconstruction  optimized  for  the  system  with  cylindrical 
transducer  array  consisting  of  rings  with  evenly  spread  piezoelectric  detectors. 

2.  METHODS 


Let  ^  to  be  a  radius  -vector  of  one  of  transducers  in  array.  The  temporal  integral  of  acoustic  pressure  detected  by 
transducer  at  the  time  t  can  be  expressed  as  the  superposition  of  waves  radiated  by  optoacoustic  sources  located  on  the 
spherical  surface  Z  of  radius  equaled  to  time  of  sound  propagation  C^t  [1]: 


u  (F,0  = 


An  Cp  c. 


-\jQ(r+cJn)dX 


(1) 


where  is  the  thermal  coefficient  of  volume  expansion  of  the  tissue;  Cp  is  the  specific  thermal  capacity.  Function 

Q{r'^  —  (r')F(r')  is  the  product  of  light  absorption  coefficient,  and  laser  fluence  F  in  the  area  where 

optoacoustic  sources  are  distributed.  Equation  ( 1 )  is  of  general  character  and  it  can  be  considered  as  the  fundamental 
imaging  equation  for  optoacoustic  tomography. 
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2.1.  Two-dimensional  solution 

In  2 -dimensional  case  equation  (1)  can  be  written  in  the  following  form: 


u{a,r)-  ^  f  ds  ,  (2) 

r{a,r) 

where  T(pc^r)  is  the  arc  of  radius  v  with  origin  in  the  center  of  receiving  transducer  with  coordinates  (coSOf,  sinOf)  . 
We  supposed  here  without  loss  of  generality  that  transducers  are  placed  on  the  cylindrical  surface  with  unity  radius.  Eq. 
(2)  presents  a  general  problem  of  integral  geometry  of  function  reconstruction  from  a  set  of  its  integrals  along  curves. 
A  common  approach  to  solution  of  this  problem  described  in  [3,4]  allows  reducing  of  eq.  (2)  to  the  set  of  simultaneous 
equations  of  Volterra  type  with  weak  singularity  kernel. 

To  define  these  equations  we  introduce  cylindrical  coordinates  x  =  p  COs6 =  siXiO  and  variable 

p  =  1  —  r  .  Let  V(OJ,  p)  =  u{(X^l  -  r)  then  following  [4]  one  can  receive  a  set  of  equations: 


^  p^+2p-p^^ 
2p 


where  (p),  (p)  are  the  Fourier  coefficients  of  corresponding  functions  v(CtjP),  /(p,0)  ■ 


•t  In  *1  -i-'t 

V*  (P)  =  —  J  e-'*“v(a,  p)da,  /,  (p)  =  —  J  f{p,d)de. 

in  ^  in  ^ 

Function  (x)  is  the  Chebyshev  polynomial  of  the  first  type.  Function  /z(p,  p)  is  expressed  as: 

(p^-2p  +  p^J 

1  + - 7 - ^  .  It  has  a  singularity  provided  that  n  — ^  n  ■ 

Ap^  -\p^  +lp- p^)  \ 


In 


Kp,p)  = 


(3) 


Figure  1.  Reduction  to  the  problem 
of  Radon  transform. 


Development  of  algorithm  of  the  spatial  distribution  of  optoacoustic 
sources  reconstruction  completely  based  on  the  exact  solution  of  equations 
(3)  has  serious  mathematical  difficulties.  In  general  we  have  to  solve  an 
infinite  set  of  equations  each  of  them  has  a  singularity.  We  have  tested 
several  methods  of  numerical  solution  of  eqs.  (3),  but  all  of  them  resulted 
in  algorithms  of  high  complexity  and  weak  convergence.  These  algorithms 
could  not  provide  a  function  reconstruction  with  necessary  precision. 

We  propose  to  use  an  iterative  approach  for  the  solution  of  the  inverse 
problem.  At  the  first  step  limited  number  of  equations  (3)  are  solved  in 
whole  area  and  this  solution  provides  an  image  with  low  resolution. 
Nevertheless  this  approximate  solution  can  be  used  for  localization  of  area 
containing  suspicious  object.  At  the  second  step  integrals  inside  area  of 
interest  are  calculated  as  difference  between  measured  data  (2)  and 
integrals  of  function  defined  at  the  first  step  over  part  of  arcs  outside  of 
area  of  interest.  At  the  third  step  we  suppose  that  area  of  interest  is  small 
compared  to  the  entire  object.  It  allows  substituting  of  the  integrals  along 
arcs  by  integrals  along  straight  lines  that  tangent  to  corresponding  arc 


curves.  As  a  result  the  problem  is  reduced  to  classical  Radon  problem.  A  filtered  backprojection  algorithm  is  employed 
for  image  reconstruction  in  the  area  of  interest  with  sufficient  resolution. 


120 


Proc.  SPIE  Vol.  4256 


2.2.  Filtered  backprojection  algorithm 

In  this  study  we  proposed  that  location  of  area  of  interest  and  integrals  of  desired  function  along  straight  lines  are 
already  known.  Therefore  ultimate  quality  of  reconstructed  images  can  be  evaluated  when  filtered  backprojection 
algorithm  is  employed. 

From  mathematical  point  of  view  the  following  integrals  along  lines  XCOS(p  +  ysin(p  =  q  defined  in  canonical 
coordinates  (^,()!))  are  known; 

+00 

u((p,q)=  I  f{qcos(p  +  ^sin(p,  qsin(p cos  (p)d^  .  (4) 

A  lot  of  approaches  are  supposed  for  the  solution  of  this  problem  [5-7].  Most  of  practical  tomography  systems  employ 
convolution  and  backprojection  algorithms  that  are  based  on  the  projection  theorem: 

w(<p,A)  =  /(Acostp,  Asincp) ,  (5) 

-t-OO 

where  ti  denotes  the  Fourier  transform  of  the  function  u  on  independent  variable  q  je^u{(p^q)  dq  . 

Two-dimensional  Fourier  transform  of  desired  function  u  calculated  in  polar  coordinates  (A^  =  A  COS(p^  A^  =  A  sin  (p)  is 

on  the  right  side  of  eq.  (5).  To  get  inversion  formula  for  the  desired  function  f  one  has  to  calculate  two-dimensional  inverse 
Fourier  transform: 


f{x,y)  =  ^  I  J  ^1  dX\  e‘^u{(p,q)dq  . 


(6) 


General  theory  of  convolution  and  backprojection  algorithms  has  been  extensively  developed  [5,8,91- 

Let  -u{(Pj^,qj)  is  the  initial  integral  defined  on  the  line  with  co ordinates .  By  definition,  coordinates 
(p^  and  qj  are  given  by  following  relations: 

(p^  =kA(p,  =  -1), 

qj  =jAq,  ;  =0,±l,±2,...iV^ 

where  Aq  is  the  discretization  step  in  radial  direction,  A(p  is  the  angular  distance  between  two  neighboring 
transducers  (see  Fig.  1).  The  cycle  condition  is  satisfied.  Value  of  Aq  can  be  evaluated  from  time 

sampling  interval  At  of  the  analog-digital  converter  employed  in  the  tomography  system:  Aq  =  At  (Cs  is  the  speed 
of  sound  in  the  tissue). 

The  phenomenological  model  of  measurement  procedure,  in  particular  the  inaccuracy  of  ultrasonic  detectors  can  be 
taken  into  account  by  introducing  apparatus  function  n(q): 

•t-oo 

-q  )u{(p,.,q')dq  ■  (s) 


A(p  = 


n 


iV, 


(p  ? 


(7) 


Algorithm  based  on  the  inversion  formula  (6)  consists  of  two  steps.  Function  Ij/ {(p ,  a)  is  computed  as  a  convolution 
of  measured  data  with  a  kernel  H^(s)  ,  and  function  /  is  resulted  as  sum  of  \f/^((p,^,  a)  over  entire  angles: 
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(9) 


J 

AL-1 


/  ^  (^,  j)  =  E  V^  ^  ft ,  ^  cos  (p,  +  y  sin  (p* )  A<p 

k=Q 

Kernel  H^(s)  depends  on  the  effective  regularizer  R‘  \A); 

H^(s)  =  ^^'fAi?^(A)cosAs  dA 
hi  ^ 


(10) 


(11) 


where  X  is  a  spatial  frequency.  A  specific  form  of  effective  regularizer  in  turn  is  defined  by  regulanzer  R{X)  and  algonthm 
of  data  interpolation  employed.  We  restricted  ourselves  to  algorithms 

described  in  the  following  way.  Let  introduce  regularizer  as  a  function,  which  satisfies  the  condition  . 

i?(0)  =  l;  i?(A)  = /?(-A);  i?(A)  =  0  if  |A|>Q, 

where  Q  =  —  is  the  Nyquist  frequency.  A  kernel  corresponding  to  the  regulanzer  is  calculated: 

Aq 

1 

(s)  =  — ^  f  A  7?  (A)  cos  Xs  dX  . 

2n  i 


It  allows  calculating  of  the  function  Wr  (<Pr  ’  at  sampling  points: 

W ni^k^nAq)  =  ^{nAq  -  jAq)u^ {(p i, ,q j)  Aq 


(12) 


Function  iu  the  intermediate  points  is  computed  by  using  a  linear  data  interpolation.  Effective  regularizer  satisfying  this 
procedure  has  a  following  form; 


sin - 


?iAq 


^  I A  +  InQ.  |i?(A  +  2nQ) 


We  used  modified  Shepp-Logan  (MSL)  7?,  (A)  and  rectangular  (RECT)  R2  (^)  spectral  windows: 


(13) 


7?,(A)  = 


sin^  ,  |A|<n;  jia)  =  P  ’ 

0,  lAl>f2; 


V  ^  7 


(14) 


Reomgula,  wmdo«.  provides  the  flat  speehal  characteristic  up  to  Nycptist  trepoency,  MSL  wmdotv  reduces  the 

”ort“st;:c™t'^ K™™  .ha,  apparams  tuuc.toh  ts  eharac.er.ced  by  uaitoru.  seusitiv.ty  aud  can  he  expressed 

in  the  form: 


niq)  = 


1  . 

0  .  W>f- 


(15) 
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3.  RESULTS 


In  this  section  results  of  the  numerical  experiments  will  be  described.  Numerical  simulations  were  performed  with 
computer  code  developed  on  the  base  of  mathematical  model  of  two-dimensional  optoacoustic  tomography  system 
(MMT).  The  basic  MMT  includes:  (1)  method  of  object  definition  ~  fanction /  ;  (2)  algorithm  of  initial  data  generation 
-  function  u  ^  ((2?^  ,q  j)  and  noise  A  ;  (3)  algorithm  of  reconstruction  -  filtered  backproj  action  and  (4)  algorithm  of  image 
quality  estimation. 

To  perform  computer  simulation  of  image  reconstruction  we  have  to  specify  the  parameters  of  the  optoacoustic 
tomography  system.  We  assume  that  a  diameter  of  the  nng  there  acoustic  transducer  are  located  is  120  mm,  number  of 

transducers  is  180,  therefore  the  angular  discretization  parameter  A(p  =2*^.  Image  reconstruction  with  the  step  of 
discretization  in  the  radial  direction  Aq=0.1  mm  is  referred  to  the  variant  1.  In  variants  2  and  3  the  calculations  were 
performed  with  steps  Aq=0.3  mm  and  Aq=0.5  mm.  Initial  data  averaged  on  the  corresponding  intervals  were  employed 
for  image  reconstruction  with  enlarged  steps: 

u^{q>„  )  =.  - - -  T  “ "  )  (rn  =  l,T).  (16) 

2m  -r  1 


We  consider  reconstructed  object  to  be  a  homogeneous  disk  of  radius  5  mm  with  characteristic  function: 


Jl  ,  (x,y)eD^-, 

jo  ,  {x,y)eDQ, 


(17) 


In  this  case  maximum  initial  signal  can  be  estimated  as  lOmm.  We  consider  the  noise  to  be  one-tenth  of  a  maximum 
signal:  o*(i/)  -  LO”^  max  (^,9)  ~  LO  mm  . 


a  b  c 

Figure  2.  Image  of  the  homogeneous  disk  of  5  mm  radius  without  noise.  The  256-grade  gray  scale  lies  in  the 
range  (-0.1  -  0.1).  a  -  Aq=0.1  mm;  b  -  Aq=0.3  mm;  c  -  Aq=0.5  mm. 


Images  of  the  homogeneous  disk  with  center  at  (-30,  0)  are  presented  in  Fig.2  and  Fig.3.  Initial  data  without  noise 
(Fig.2)  and  with  noise  (Fig.3)  were  used  for  the  image  reconstruction.  Calculations  were  performed  for  three  variants 
corresponding  to  different  steps  of  discretization  in  radial  direction  Aq  in  the  square  with  60-mm  side.  MSL  spectrum 
window  was  employed. 

The  tomograms  presented  in  Fig.  2  allows  to  estimate  the  artifacts  level  for  tree  variants  of  reconstruction.  By  artifact 
one  usually  mean  different  kinds  of  systematic  reconstruction  errors.  Unlike  noise,  artifacts  demonstrate  a  clear  well- 
ordered  geometric  structure  that  depends  on  the  nature  of  systematic  errors.  The  artifacts  presented  in  Fig.  2  are  the 
result  of  discretization  errors.  Their  ray-like  structure  is  defined  by  the  significant  difference  in  the  numbers  of  angular 
and  spatial  samples.  Artifact  level  can  be  estimated  as  a  ratio  of  maximum  deviation  on  the  tomogram  outside  of  the 
disk  image  to  the  mean  value  inside  disk.  We  received  the  following  values:  variant  VI  (Aq=0.1  mm)  -  10%,  V2 
(Aq=0.3  mm)  -  4.4%,  V3  (Aq=0.5  mm)  -  2.7%. 
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In  practical  tomography  the  discretization  artifact  levels  are  small  and  dominated  by  noise  and  other  errors.  It  is  clearly 
demonstrated  in  Fig.3  where  the  image  of  the  same  object  in  the  presence  of  noise  is  shown.  Noise  level  is  three  times 
higher  than  artifacts,  but  in  somewhere  the  regular  structure  is  slightly  appeared.  Noise  level  can  be  defined  as  mean 
square  deviation  on  the  tomogram  calculated  in  the  region  sufficiently  far  from  the  object.  The  noise  levels  determined 
for  three  variants  of  image  reconstruction  are  presented  in  Table! : 

Table  1 .  Parameters  of  the  filters. 


VI 

V2 

V3 

MSL 

0.304 

0.058 

0.027 

RECT 

0.657 

0.126 

0.059 

a  b  c 

Figure  3.  Image  of  the  homogeneous  disk  of  5  mm  radius  with  noise.  The  256-grade  gray  scale  lies  in  the  range  (-0.3 

-  03).  a  -  Aq=0.1  mm;  b  -  Aq=03  mm;  c  -  Aq=0.5  mm. 


Noise  level  is  significantly  reduced  then  averaging  is  applied,  the  tomogram  becomes  smoother.  But  in  the  same  time 
sharpness  of  the  object  edges  is  also  reduced,  and  the  spatial  resolution  decreases.  To  demonstrate  the  special  resolution 
provided  by  the  filtered  backprojection  algorithm  the  images  of  two  disks  separated  by  a  small  gap  were  reconstructed. 
Calculations  were  performed  with  noisy  initial  data  and  MSL  spectral  window  was  employed.  Results  are  shown  in 
Fig. 4.  The  0.1  mm  gap  is  clearly  resolved  only  on  the  tomogram  calculated  with  step  Aq=0.1  but  diis  tomogram  is 
noisy  one.  Calculations  with  variants  2  and  3  provide  more  smooth  tomogram  but  with  low  resolution.  The  03  mm  gap  is 
resolved  on  entire  tomograms. 


4.  DISCUSSION 

The  main  problem  discussed  in  this  section  is  the  choice  of  quantitative  parameters  describing  the  quality  of 
homographic  images.  In  practice,  this  problem  is  solved  by  choosing  several  basic  quality  criteria  and  describing  the 
corresponding  measurement  methods.  We  will  use  following  quality  criteria:  (1)  noise  level  on  the  tomogram,  (2)  jump 
transfer  function,  (3)  loss  contrast  function  and  (4)  contrast-dimension  curve.  These  criteria  allow  us  to  optimize  the 
reconstruction  algorithm  and  compare  the  quahty^  of  images  obtained  in  various  tomography  systems.  General  theory 
developed  in  [9]  provides  theoretical  estimations  of  the  main  parameters  of  reconstruction  quaiity^  We  will  adapt  them 
for  opto  acoustic  tomography  system  and  compare  theoretical  estimations  with  results  of  numerical  simulation. 

Noise  level  on  the  tomogram  can  be  estimated  as 
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(18) 


, ,  ^  a"(w)A<p 

o-(/'  )  =  • 

An 


dA 


where  g{u)  is  the  mean  square  deviation  of  die  initial  data  noise  and  (  a)  is  the  effective  reguianzer.  For  the  algorithm 
we  used  in  our  calculation  this  formula  can  be  simplified: 


Aq=0.3  mm 


;Mii 


a  b 

Figure  4.  Images  of  two  disks  of  radii  5  mm  and  3  mm  with  a  small  gap  between  them.  The  gap  value  is  0.  i  mm  (a) 
and  0.3  mm  (b).  The  256-grade  gray  scale  lies  in  the  range  (-0.5  -  1 .5). 


^q- 


(19) 


where  constant  value  {A)  depends  on  spectrum  window: 
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(20) 


_f  1/1271==  0.02658  ,  R  =  R,; 

^^*''^^~|7r/18-l/67r  =  0.121481  ,  R  =  R^. 

In  particular,  for  the  case  when  steps  of  discretization  Acp  =  2^  and  Aq=0.1  mm  (variant  VI),  spectral  window  R1  and 
initial  data  with  noise  (j(w)  =  0.1  cm  ,  formula  (15)  give  the  30%-noise  on  the  tomogram.  This  value  corresponds  well  to 
results  of  our  numerical  simulations  presented  in  Table  1. 

A  jump  transfer  function  (JTF)  can  be  used  as  a  measure  of  spatial  resolution  on  the  resulted  image.  The  JTF 
characterizes  the  smoothing  of  discontinuity  of  object  edges  on  the  tomogram.  To  define  JTF  one  has  to  restore  the 

characteristic  function  of  the  disk  with  radius  b  »  Aq  provided  /"^  =  1  in  center  part  of  the  disk.  Let  the  disk  center 
has  coordinates  ,  then  reconstruction  of  function  {nAx^Uy)  with  the  small  step  Ax«Aq  results  in  JTF: 

(x)  =  (x  -f  <3^  H-  Z?,  )  X  >  0  .  Theoretical  estimation  provides  the  following  formula  for  JTF: 

p^(jc)  =  (21) 

2  7t{  A 

where  n(A)  is  the  Fourier  transform  of  apparatus  function. 

Figure  5  shows  jump  transfer  functions  corresponding  to  MSL  spectral  window  (three  upper  curves)  and  RECT  window 


x/Aq 


(three  lower  curves).  These  data  were  resulted  from  reconstruction  of  the  characteristic  function  of  the  disk  with  radius 
5  mm  without  noise.  Calculations  were  performed  along  X-axis  in  the  range  [-25.0,  -25.0+2Aqi]  with  step 

Ax  ^  Aq.  /  200  (i=l,2,3  correspond  to  variant  of  reconstruction).  Rectangular  window  provides  more  sharp  edge  of 
image;  at  x=0  as  it  follows  from  formula  (16)  these  curses  have  the  same  value  equaled  to  0.5. 
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Contrast  loss  function  can  be  evaluated  from  the  reconstruction  of  the  small  disks.  Fig.6  represents  tomogram  profiles 
along  x-axis  obtained  in  exactly  the  same  way  as  a  jump  transfer  function,  but  with  disks  of  different  radii:  b=2Aq=0.2 
mm,  b=Aq,  b=Aq/2  and  with  the  x-variable  representing  the  distance  from  the  center  of  the  disk.  Disk  center  was  at  the 
point  (-30.0,  0). 


c 


0.8 


0.4 


0.0 


Figure  6.  Images  of  small  disks.  Curves  correspond  to  the  following  parameters:  1  -  b=0.2  mm,  VI,  RECT;  2  - 
b=0.2  mm,  VI,  MSL;  3  -  b=0.1  mm,  V2,  RECT;  4  -  b=0.1  mm,  V2,  MSL;  5  -  b=0.05  mm,  V3,  RECT;  6  - 
b=0.05  mm,  V3,MSL. 

The  even  curves  were  reconstructed  with  MSL  window;  they  are  flatter  than  corresponding  curves  reconstructed  with 
RECT  window.  Contrast  C  in  the  disk  image  center  is  reduced  with  radius  decreasing.  We  used  the  contrast  definition 

as  excess  of  function  over  background  related  to  the  background  itself:  C  ~  — —I  -  Dependence  of 

a  G 

contrast  in  the  center  of  a  small  object  on  its  dimension  is  described  by  the  contrast  loss  function  (b)  .  Theory 
provides  the  following  estimation  for  the  contrast  loss  function  in  the  case  of  small  disk  with  radius  b  [9]: 


oo 

C^{b)^b\n{X)R^{X)J,{Xb)dX,  (22) 

0 

where  is  Bessel  function  of  the  1-st  order.  Fig. 7  shows  the  contrast  loss  function  calculated  for  variants  V1-V3 

with  MSL  window. 

The  contrast-dimension  function  {b)  is  defined  from  the  relation: 

CD\b)^C\b)^Ycy{f")  (23) 

where  G  is  the  noise  level  on  the  tomogram,  y  is  a  given  number  called  as  the  detection  level.  As  follows  from  this 
definition  the  {b)  is  the  contrast  value  such  that  the  corresponding  signal  is  Y  times  more  than  the  noise  level. 


Proc.  SPIE  Vol.  4256 


127 


Figure  7. Contrast  loss  function. 


Usually  the  detection  level  is  chosen  as  3.  Contrast-dimension  curve  divides  the  plot  surface  on  two  parts;  area  of 
(upper  the  curve)  .»d  iuvieible  cues  (below  the  curve).  Oubast-dimeusiou  funcou  ts  a  un^e^al 
characteristic  of  any  tomography  system.  It  can  be  employed  for  X-ray,  MRI,  Ultrasound  systems  for  their  performance 
co™S  fi^Tebows  JiudLon  of  0.2  mm  radius  can  be  revealed  on  the  tomogram  provded  ,ts  contrast  ,s  1  o, 
higher.  For  smaller  objects  (0.1  mm,  for  example)  the  contrast  should  be  3  and  higher. 


Disk  radius,  mm 

Figure  8.  Contrast-dimension  function  when  rectangular  spectrum  window  was  used. 


It  IS  necessary  to  emphasize  that  the  functions  {x\C\b)  charactenzing  the  spatial  resolution  depends  on  only  ratios 
x/Ao,  i?/Aa-  Therefore  the  expressions  (17),  (18)  derived  for  these  functions  have  a  universal  character,  if  all 
dimensions  are  normalized  by  discretization  step  Aq.  It  allows  us  to  employ  results  of  [9]  and  this  work  for  optoacoustic 
systems  of  various  configurations. 
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5.  CONCLUSION 


This  work  is  the  first  step  for  a  development  of  mathematical  model  of  image  reconstruction  for  a  two-dimensional 
optoacoustic  imaging  system  with  cylindrical  transducer  array.  We  considered  only  one  stage  of  this  model  namely  filtered 
backprojection  algorithm  in  assumption  that  integrals  along  straight  lines  are  known.  We  are  planning  to  generalize  this 
algorithm  to  the  case  of  integrals  along  arcs.  It  will  allow  us  to  reconstruct  images  of  the  inclusions  located  in  the  transducers 
vicinity. 
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ABSTRACT 

Thermal  mechanism  is  fundamental  for  the  generation  of  ultrasonic  waves  in  a  course  of  absorption  of  laser  radiation.  Its 
efficiency  is  relatively  low,  but  increases  significantly  with  temperature  (thermal  non-linearity).  If  the  phase  transition  of  the 
irradiated  medium  occur,  the  efficiency  may  exceed  its  linear  regime  value  several  orders  oT  magnitude.  The  idea  of 
optoacoustic  supercontrast  for  early  cancer  detection  is  based  on  this  fact.  We  performed  feasibility  studies  and  describe 
requirements  to  and  properties  of  the  optoacoustic  supercontrast  agent  based  on  nanoscopic  particles.  The  results  of  the 
preliminary  experiments  with  the  metal  and  carbon  nanoparticles  as  optoacoustic  supercontrast  are  presented.  Theoretical 
estimations  show  the  possibility  to  enhance  the  efficiency  of  optoacoustic  generation  up  to  three  orders  of  magnitude  under 
irradiation  conditions  of  laser  optoacoustic  imaging  in  the  depth  of  human  tissue. 

Keywords:  optoacoustic  effect,  light  absorption  and  scattering,  optoacoustic  transients,  nanoparticles 


1.  INTRODUCTION 

The  optoacoustic  tomography  (OAT)  is  a  novel  promising  technique  for  diagnostic  imaging  of  cancer  [1-4].  OAT 
employs  illumination  of  tissue  under  study  with  short  laser  pulses,  generation  of  ultrasonic  waves  in  tissue  due  to 
absorption  of  laser  radiation  and  detection  of  profiles  of  ultrasonic  waves  with  high  temporal  resolution  by  an  array  of 
ultrawide-band  acoustic  transducers  [5,6].  The  optoacoustic  imaging  takes  advantage  of  high  optical  contrast  between 
normal  and  cancerous  tissue  and  low  distortion  of  ultrasonic  waves  delivering  information  on  the  distribution  of 
absorbed  optical  energy.  Estimations  show,  that  it  is  possible  to  detect  an  absorbing  sphere  with  I  mm  in  diameter  and 
optical  properties  of  cancerous  tissue  at  a  depth  of  7  cm  in  the  breast  [7,8].  However,  the  exact  value  of  optoacoustic 
contrast  in  small  early  tumors  is  yet  to  be  determined.  Currently,  it  is  believed  that  the  main  portion  of  the  optoacoustic 
contrast  comes  from  tumor  angiogenesis.  It  is  unlikely  that  an  advanced  angiogenesis  can  be  developed  at  the 
microscopic  stage  of  early  tumors.  Therefore,  further  enhancement  of  optoacoustic  imaging  contrast  will  be  useful  for 
detection  of  early  cancer  at  significant  depth.  Simultaneously,  the  enhanced  contrast  would  permit  application  of  much 
lower  laser  energy  for  irradiation,  which  in  turn  may  replace  flash-lamp  pumped  lasers  currently  employed  in  OAT 
with  inexpensive  and  compact  diode  lasers.  In  addition,  greater  sensitivity  of  OAT  will  allow  the  use  of  greater 
repetition  rate  of  laser  pulses,  yielding  greater  signal-to-noise  ratio.  The  optoacoustic  contrast  is  determined  not  only 
by  the  variations  of  light  absorption  and  scattering  coefficient,  but  also  by  the  variation  of  the  optoacoustic  efficiency. 
Thus,  the  idea  is  not  to  enhance  average  light  absorption  in  tissue  (as  conventional  exogenous  contrast  agents  would), 
but  to  enhance  the  optoacoustic  efficiency.  It  is  well  established,  that  generation  of  acoustic  waves  is  substantially 
more  effective  if  the  evaporation  or  other  gas  production  takes  place  [9-13],  On  the  other  hand,  for  the  optoacoustic 
imaging  of  biological  tissue  to  be  safe,  the  average  temperature  rise  may  not  exceed  1°C.  Therefore,  the  evaporation 
can  take  place  only  at  extremely  small  absorbing  centers  in  the  tissue  not  to  damage  the  tissue.  To  achieve  threshold 
for  evaporation,  photo  dissociation  or  other  gas  production  at  a  safe  level  of  laser  energy  fluence,  the  absorbing  particles 
should  have  dimensions  on  the  nanometer  scale.  The  same  limitation  for  the  particle  size  comes  with  necessity  for  the 
contrast  agent  to  diffuse  between  ceils,  and  penetrate  cell  membranes. 


2,  THEORY  OF  OPTOACOUSTIC  PHENOMENA  IN  HETEROGENEOUS  MEDIUM 

Leu  us  for  simplicity  reason  consider  a  turbid  medium  as  a  clear  medium  with  immersed  absorbing  particles  (R^  -  radius  of 
the  particle).  Let  concentration  of  the  particles,  n,.  be  sufficiently  high,  so  that  the  quantity  of  the  particles  in  the 
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elementary”  volume  of  a  x  a  x  a  (see  Fig.l)  is  substantially  greater  then  unity:  nc^  »  1. 


Figure  1.  Schematics  of  heterogeneous  medium  with  distributed  nanoparticies. 

On  the  other  hand,  the  concentration  should  be  limited,  so  that  the  optical  absorption  over  the  “elementary”  volume  can  be 

taken  as  homogeneous:  (yna  «  1  =  -  cross-section  of  optical  absorption  of  the  particle).  Particle  absorbs  laser 

radiation  and  it  temperature  rise  up.  Surrounding  liquid  is  heated  due  to  heat  diffusion  from  the  particle.  Pulsed  laser  heating 
of  absorbing  particle  immersed  in  transparent  liquid  can  be  treated  in  several  stages. 

If  the  temperature  of  the  particle  does  not  exceed  the  boiling  point  of  the  liquid,  (lower  limit  of  the  laser  fluence),  the 
phase  transition  cannot  occur. 


Figure  2,  The  system  of  coordinates  for  calculation  of  a  small  volume  irradiation. 

The  temperature  of  the  liquid  around  the  particle  after  irradiation  by  the  laser  pulse  can  be  taken  in  the  following  form  (the 
temperature  fields  of  the  particles  don’t  overlap  due  to  short  laser  pulse  duration  and  low  concentration  of  the  particles): 

T  =  TqHTp-  ro)-^exp(-^~-^^^ );  r>i?o.  (1) 

where  X  is  diffusivity  of  the  liquid,  is  laser  pulse  duration,  Tp  <  Tg^  is  the  temperature  of  the  particle  (we  consider  the 
temperature  in  the  particle  homogeneous). 
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According  to  (1),  the  temperature  of  the  particles  can  be  found  as: 


T  =T  + 


4(Pc^),-7^(V^  +  2j:  +  C/x)’ 


(2) 


where  p  is  the  density  of  the  medium  (subscript  L  stands  for  “liquid”,  and  subscript  p  stands  for  “particle”),  is  the 


specific  heat,  x  =  X'^L  ^  ^  is  the  ratio  of  heat  diffusion  length  and  particle  radius,  ^  = 
heat  capacity  of  the  particle  and  surrounding  liquid. 


_  (PCp)p 


XPCp)l 


characterizes  the  ratio  of 


Evaporation  starts  when  the  temperature  of  the  particle  reaches  the  boiling  point.  This  threshold  of  evaporation  can  be 
estimated  as: 

The  threshold  has  a  minimum  at  certain  particle  radius,  ^min  ~  ^  2  .  This  gives  an  optimal  radius  of  the  particle  In  order 


to  minimize  the  evaporation  threshold: 


^opt 


I6(pcp)i 

j  ipCp)p 


X^L  ■ 


(4) 


The  optimal  radius  of  the  particle  depends  on  the  laser  pulse  duration,  specific  heat  of  the  particle  and  and  specific  heat  of  the 
surrounding  liquid.  For  carbon  or  metal  particle  in  water  and  10  ns  laser  pulse  optimal  radius  is  of  the  order  of  110  nm  to  130 
nm.  As  thermal  diffusion  length  in  the  particle  is  of  the  order  of  1  /xm,  the  model  of  homogeneously  heated  particle  is  valid. 
In  these  conditions  the  threshold  of  evaporation  may  be  presented  as: 

In  an  optimum  situation,  heat  containing  in  the  surrounding  liquid  4-5  times  exceeds  that  of  the  particle.  If  the  laser  fluence 
exceeds  the  threshold  (5),  evaporation  at  the  particle  takes  place.  If  laser  fluence  does  not  exceed  the  threshold  (5),  then  the 
evaporation  of  the  liquid  does  not  occur  and  excitation  of  ultrasound  wave  takes  place  due  to  thermal  expansion  of  the 
particle  and  surrounding  layer  of  liquid.  As  the  local  temperature  rise  can  be  several  tens  of  degrees  and  higher,  the  so-called 
thermal  nonlinearity  phenomema  can  occur.  The  thermal  nonlinearity  leads  to  the  rise  of  the  efficiency  of  optoacoustic 
generation  3-4  times  in  a  case  of  water  (see  below). 

If  the  laser  fluence  exceeds  threshold  (5)  the  temperature  of  the  particle  and  surrounding  liquids  comes  to  boiling  point  at 
some  moment  in  a  course  of  laser  pulse.  Then  thin  vapor  layer  is  generated  around  the  particle.  Due  to  lower  heat 
conductivity  of  vapor  (relatively  to  that  of  Liquid),  the  temperature  of  the  particle  continues  to  increase.  After  the  laser  pulse 
the  temperature  of  the  particle  relax,  to  the  initial  value  due  to  heat  fluence  from  the  particle  through  the  vapor  layer  into  the 
liquid.  The  total  mass  of  the  vapor  is  determined  by  the.  heat  transmitted  from  the  particle  to  the  liquid.  In  a  course  of 
relaxation  at  some  moment  of  time  the  temperature  of  the  particle  becomes  lower  then  the  boiling  point  and  condensation  of 
the  vapor  starts.  Total  mass  of  vapor  produced  around  single  particle  does  not  exceed  the  foOowing  value: 

(6) 

where  X  is  the  heat  of  evaporation.  We  neglect  the  heat  contain  of  the  particle. 

The  pressure  wave,  generated  by  heated  elementary  volume  can  be  taken  in  the  form: 

(7) 

A7tr  dt^ 

where  (51/  is  a  variation  of  the  elementary  volume  due  thermal  expansion  and  vapor  production. 
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In  case  of  homogeneous  bulk  absorption  of  light  the  expansion  of  an  elementary  volume,  5V ,  can  be  expressed  as  in  [7]: 


/3 


0 


(8) 

a-'  pcp 

where  |3  is  the  volume  expansion  coefficient.  Thermal  expansion  coefficient  depends  on  the  temperature  and  for  water  this 
dependence  is  quite  strong  (so-called  “thermal  nonlinearity”,  see  Fig.3).  Therefore,  before  the  temperature  of  the  particle  can 
reach  the  boiling  point,  the  efficiency  of  the  optoacoustic  generation  increases  up  to  2-3  times. 


Figure  3.  Thermal  expansion  coefficient  of  water  as  a  function  of  temperature. 
The  thermal  sources,  producing  the  acoustic  wave  can,  taken  as  follows: 

dp 


8T 


W- 


(9) 


Expression  (9)  shows,  that  for  low  temperature  rise  j/^Q  ^  d^jdT  5t|  « 1  the  first  term  in  the  brackets  dominates  and  the 
distribution  of  heat  source  is  proportional  to  the  temperature  rise  ST .  In  an  opposite  case  of  relatively  high  temperature  rise 
j^Q  ^  dP  jdT  ST  »\  the  second  term  is  much  greater  then  the  first  and  the  distribution  of  heat  source  is  proportional  to 

\  Thus  for  short  laser  pulse  the  spatial  distribution  of  heat  sources  will  be  a  sum  of  two  exponents: 


PCp 


pc 


P  ) 


^^-IpLaZ 

dT 


(10) 


For  low  temperature  rise  (low  fluence)  the  leading  edge  of  the  optoacoustic  signal  (that  resembles  the  spatial  distribution  of 
heatisources)  will  have  exponential  shape  with  the  exponent  higher  fluence  the  exponential  slope  will  be  2  times 

sharpen  2/J.QCQt ,  In  a  case  of  heterogeneous  medium  with  absorbing  centers,  the  variation  of  the  elementary  volume  due  to 
vapor  production  can  be  estimated  as  (we  consider  pressure  relaxed) : 


^2  ^  l^a^O 

PvA 

where  is  the  vapor  density.  Heating  of  the  vapor  will  give  further  rise  to  the  elementary  volume: 

sv^  _  r-1 

PO  y 


(11) 


(12) 
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where  y  is  the  adiabatic  exponent,  pQ  is  ambient  pressure. 


Let’s  compare  the  efficiency  of  various  mechanisms  of  the  optoacoustic  phenomena.  Volume  deformation  per  unit  specific 
heat  release  for  various  mechanisms  are  presented  in  the  Table  below. 


Thermal  expansion 
of  liquid  water 

Vapor 

Production 

Thermal  expansion 
of  water  vapor 

SVj ,  crn^  IJ 

8.6*10"^  n 

3.2*10"^ 

1.7 

Thus,  light  absorption  by  submicron-size  particles  can  enhance  the  efficiency  of  optoacoustic  generation  by  several  orders  of 
magnitude,  but  this  process  has  a  threshold  dependent  on  laser  pulse  duration  and  size  of  the  particle. 


3.  MATERIALS  AND  EXPERIMENTAL  SET-UP 

Studies  of  the  optoacoustic  phenomena  in  heterogeneous  media  were  made  with  the  Q-switched  Nd-YAG  laser  (laser  pulse 
duration  ^14  ns  at  1/e  level)  operated  at  fundamental  ( =  1 .064  fXm )  wavelength  and  the  second  harmonic 
( =  0.532  \X7n).  The  laser  pulse  energy  was  varied  by  neutral  density  filters.  Pulse  repetition  rate  did  not  exceed  3  Hz  and 
was  chosen  to  fulfill  the  safety  requirement  not  to  elevate  the  average  temperature  of  the  medium  more  than  1  degree. 
Actually,  most  of  the  measurements  have  been  done  with  single  laser  pulses. 

We  used  optically  heterogeneous  media:  water  solution  of  India  ink  of  various  concentrations,  water  suspension  of  gold 
nanoparticles  (diameter  10  nm,  30  nm  and  100  nm)  with  the  volume  concentration  of  the  particles  of  the  order  of  10  .  As  a 
reference  homogeneous  medium  water  solution  of  cupric  chloride  of  various  concentrations  was  employed. 

The  enhancement  of  the  optoacoustic  imaging  contrast  in  heterogeneous  medium  was  investigated  with  laser  optoacoustic 
imaging  system  (LOIS-2)  designed  for  breast  cancer  imaging.  The  system  uses  32-element  arc-shaped  array  of  piezoelectric 
transducers  with  ultrasonic  frequency  band  of  0.05-3  MHz,  which  provides  spatial  resolution  of  the  order  of  0,5  mm. 


4.  EXPERIMENTAL  RESULTS 

The  leading  edges  of  optoacoustic  signals  generated  in  1%  solution  of  India  ink  in  water  for  various  energy  fluences  are 
depicted  in  Fig.  4.  The  light  absorption  coefficient  of  this  solution  was  /i^  =  18  cm  . 


Figure  4.  Profiles  of  leading  edges  of  the  optoacoustic  signals  in  aqueous  solutions  of  1%  India  ink  generated  by  laser  pulses 
with  three  different  fluences. 
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For  low  fluence  (<5*o  ~  ^ I cm^)  the  leading  edge  has  exponential  shape  with  the  exponent,  consistent  with  the  light 

absorption  coefficient  (see  (10)).  With  an  increase  of  the  incident  fluence  (Oq  =  88.7  mJ/ cm  ),  the  slope  becomes  steeper 

and  corresponds  to  2  times  sharper  exponent  (see  (10),  the  second  term).  In  the  intermediate  case  {^q  =32 A  mJ I cm^) 
beneath  the  surface  the  slope  is  two  times  higher,  then  in  the  depth  of  the  irradiated  volume  (see  Fig.4).  This  is  in  consistence 
with  the  “thermal  nonlinearity”  model  (see  (9)-(10)). 

Let’s  discuss  the  intermediate  case  Oq  =  32.4  mJ/cm^  in  more  details.  If  we  assume  homogeneous  absorption  of  light  in 
the  aqueous  solution,  the  maximum  temperature-rise  will  be  AT  =  pc ^  =0.14K.  On  the  other  hand,  the  influence 

of  thermal  nonlinearity  is  so  significant,  that  it  canot  be  expected  for  such  a  low  temperature  elevation.  The  local 
temperature-rise  can  be  estimated  from  the  ratio  of  the  amplitudes  of  exponents  in  (10).  Fit  of  the  leading  edge  of  the 

optoacoustic  signal  with  the  formula  (10)  yields  -^ST  =  0.51,  which  makes  it  possible  to  estimate  the  local 

temperature  elevation  as  ST  =  10  .K”  (for  water  at  room  temperature  d^jd-T  =  0.05  K  ^ ).  Thus,  the  absorbing  carbon 
particles  are  exposed  to  at  least  70  times  higher  temperature-rise  then  the  average  temperature  elevation  in  the  irradiated 
volume  of  the  solution.  Inside  the  practice  the  temperature  is  even  higher. 


Figure  5.  Amplitude  of  optoacoustic  signal  vs.  energy  fluence  for  water  suspension  of  carbon  particles. 

The  dependence  of  the  amplitude  of  optoacoustic  signal  vs.  laser  energy  fluence  for  water  suspension  of  carbon  particles  is 
presented  in  Fig.5.  The  diameter  of  the  particles  was  close  to  optimal  (see  (4)),  so  it  was  reasonable  to  expect  mamfestation 

of  evaporation  at  absorbing  centers  for  the  laser  fluence  of  approximately  (15  -  20)  and  significant  thermal 

nonlinearity  below  this  limit. 

The  presented  data  is  consistent  with  the  theoretical  model  discussed  above.  For  fluences  ^Q<SmJlcm  the  amplitude  of 
the  optoacoustic  signal  generated  in  heterogeneous  medium  (points  in  Fig.5)  is  proportional  to  the  fluence  and  coincides  well 
with  the  theoretical  value  for  homogeneous  medium  (solid  line).  This  means  that  the  temperature-rise  at  the  absor  mg 
particle  is  relatively  low  (at  least  below  10  K)  and  no  thermal  nonlinearity  is  manifested. 

In  the  fluence  range  of  S  mJ  /  cm^  <^o  <15  mJ /  cm^  the  amplitude  of  the  optoacoustic  signal  generated  in  heterogeneous 
medium  exceeds  that  in  a  homogeneous  medium  up  to  3  times.  This  corresponds  well  to  the  rise  of  optoacoustic  generation 

efficiency  due  to  thermal  nonlinearity  below  the  evaporation  threshold.  For  higher  fluence  q  >  15  mJ  /  cm  a  sharp 
increase  of  the  amplitude  of  excited  signals  takes  place.  It  may  occur  due  to  the  elevation  of  the  temperature  of  particles  to 
the  boiling  point  and  following  evaporation  of  water  layer  around  the  absorbing  particle.  In  accordance  with  the  enhancement 
of  the  efficiency  of  optoacoustic  generation  in  vapor  production  regime  (see  Table  1),  a  dramatic  increase  of  the  amplitude  o 
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Optoacoustic  images  of  the  same  vessels  filled  with  the  water  suspension  of  gold  particles  of  various  diameter  at 
A  — 0.53umare  presented  in  Fig. 7  (volume  concentration  of  the  particles  was  3*10'^^.  For  the  fluence  below  the 
evaporation  threshold  (Fig. 7a)  the  100  nm  in  diameter  particles  i  vessel  J)  produce  the  brightest  image,  due  to  higher  light 
absorption  coefficient  of  the  solution.  Correspondingly,  the  solutions  of  30  nm  (vessel  2)  and  10  nm  (vessel  1)  particles 
produce  the  images  with  lower  intensity.  For  the  fluence  above  evaporation  threshold  (Fig. 7b)  the  brightest  image  was 
produced  by  10  nm  particles,  and  the  intensity  of  the  images  is  reverse  to  that  below  the  threshold. 


5.  DISCUSSION  AIND  CONCLUSIONS 

Light  absorbing  nanoparticles  were  studied  as  a  potentially  very  effective  contrast  agent  for  laser  optoacoustic  imaging  in 
tissue.  Enhanced  efficiency  of  acoustic  wave  generation  can  take  place  due  to  superheating  of  absorbing  particles  followed 
by  evaporation  of  the  surrounding  water  layer.  Some  other  phenomena  of  phase  transition  and  gas  production  can  be  utilized. 
Laser  generation  of  ultrasound  in  heterogeneous  medium  can  be  2-3  orders  more  efficient  than  that  in  optically  homogeneous 
medium.  An  optimal  particle  size  yielding  minimal  threshold  of  evaporation  at  the  particle  surface  was  found  in  the  range  of 
-'100-nm.  The  threshold  of  evaporation  is  relatively  high  and  exceeds  10  mJ/cm^  for  metal  or  carbon  particles  suspension  in 
water.  However,  application  of  the  plasmon  resonance  absorption  in  the  particle,  the  threshold  can  be  lowered  drastically 
(due  to  higher  crossection  of  light  absorption).  One  can  expect  a  decrease  of  the  threshold  up  to  2-3  orders  of  magnitude  in 
case  of  utilization  of  plasmon  resonance  for  nanoparticie  absorption. 

Nanotechnology  is  a  rapidly  expanding  field  of  science  and  engineenng  that  offers  a  unique  potential  for  integrating  contrast 
agents  and  receptor-specific  markers  for  cancer,  promising  to  yield  smart  contrast  agents.  Gold  nanoshells  filled  with  silica 
can  be  produced  to  absorb  strongly  the  desired  wavelength  of  near  infrared  laser  pulses  [14],  These  nanoparticles  may  be 
further  modified  to  yield  exponential  enhancement  of  optoacoustic  signals  upon  pulsed  laser  irradiation.  Effective 
bioconjugates  have  been  found  to  link  various  nanoparticies  with  biological  molecules  [15].  The  problem  of  selective 
targeting  of  cancerous  ceils  and  tissues  is  currently  also  under  active  development.  Peptides,  other  ligands  and  antibodies 
specific  to  surface  receptors  in  breast  cancer  cells  or  endothelial  cell  of  tumor  microcapillary  network  are  available  for 
selective  targeting  of  malignant  tumors  [16,17].  Thus,  all  components  are  available  for  successful  development  of  the 
optoacoustic  imaging  technique  with  supercontrast.  A  receptor- specific  contrast  agent  in  combination  with  optoacoustic 
imaging  may  dramatically  enhance  the  level  of  sensitivity  of  detection  of  small  breast  tumors,  reduce  the  rate  of  false 
positives  and  negatives,  provide  diagnostic  imaging  information,  and  be  useful  in  monitoring  for  early  stages  of  breast  cancer 
and  cancer  recurrence. 

We  conclude  that  absorbing  nanoparticies  conjugated  with  targeting  agents  specific  to  cancer  may  be  applied  as  a  particulate 
contrast  agent  for  early  detection  of  cancer  with  laser  optoacoustic  imaging  system. 
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ABSTRACT 


Laser  induced  phase  transition  in  the  liquid  surrounding  absorbing  nanoparticle  results  in  dramatic  enhancement  of  the 
thermoacoustic  efficiency  and  excitation  of  high-amplitude  pressure  waves.  Feasibility  study  were  performed  in  tissue 
phantoms.  Experiments  supported  with  theoretical  model  demonstrated  possibility  to  increase  the  efficiency  of  optoacoustic 
generation  in  tissue  up  to  three  orders  of  magnitude  by  application  of  absorbing  nanoparticles  with  optimal  diameter. 

Keywords:  optoacoustic  effect,  light  absorption  and  scaitering,  optoacoustic  transients,  nanoparticles 


1.  INTRODUCTION 

The  optoacoustic  tomography  (OAT)  is  a  novel  promising  technique  for  diagnostic  imaging  of  cancer  [1-4].  OAT 
employs  illumination  of  tissue  under  study  with  short  laser  pulses,  generation  of  ultrasonic  waves  in  tissue  due  to 
absorption  of  laser  radiation  and  detection  of  profiles  of  ultrasonic  waves  with  high  temporal  resolution  by  an  array  of 
ultrawide-band  acoustic  transducers  [5,6].  The  optoacoustic  imaging  takes  advantage  of  high  optical  contrast  between 
normal  and  cancerous  tissue  and  low  distortion  of  ultrasonic  waves  delivering  information  on  the  distribution  of 
absorbed  optical  energy.  Estimations  show,  that  it  is  possible  to  detect  an  absorbing  sphere  with  1  mm  in  diameter  and 
optical  properties  of  cancerous  tissue  at  a  depth  of  7  cm  in  the  breast  [7,8].  However,  the  exact  value  of  optoacoustic 
contrast  in  small  early  tumors  is  yet  to  be  determined.  Currently,  it  is  believed  that  the  main  portion  of  the  optoacoustic 
contrast  comes  from  tumor  angiogenesis.  It  is  unlikely  that  an  advanced  angiogenesis  can  be  developed  at  the 
microscopic  stage  of  early  tumors.  Therefore,  further  enhancement  of  optoacoustic  imaging  contrast  will  be  useful  for 
detection  of  early  cancer  at  significant  depth.  Simultaneously,  the  enhanced  contrast  would  permit  application  of  much 
lower  laser  energy  for  irradiation,  which  in  turn  may  replace  flash- lamp  pumped  lasers  currently  employed  in  OAT 
with  inexpensive  and  compact  diode  lasers.  In  addition,  greater  sensitivity  of  OAT  will  allow  the  use  of  greater 
repetition  rate  of  laser  pulses,  yielding  greater  signal-to-noise  ratio.  The  optoacoustic  contrast  is  determined  not  only 
by  the  variations  of  light  absorption  and  scattering  coefficient,  but  also  by  the  variation  of  the  optoacoustic  efficiency. 
Thus,  the  idea  is  not  to  enhance  average  light  absorption  in  tissue  (as  conventional  exogenous  contrast  agents  would), 
but  to  enhance  the  optoacoustic  efficiency.  It  is  well  established,  that  generation  of  acoustic  waves  is  substantially 
more  effective  if  the  evaporation  or  other  gas  production  takes  place  [9-13],  On  the  other  hand,  for  the  optoacoustic 
imaging  of  biological  tissue  to  be  safe,  the  average  temperature  rise  may  not  exceed  1°C.  Therefore,  the  evaporation 
can  take  place  only  at  extremely  small  absorbing  centers  in  the  tissue  not  to  damage  the  tissue.  To  achieve  threshold 
for  evaporation,  photodissociation  or  other  gas  production  at  a  safe  level  of  laser  energy  fluence,  the  absorbing  particles 
should  have  dimensions  on  the  nanometer  scale.  The  same  limitation  for  the  particle  size  comes  with  necessity  for  the 
contrast  agent  to  diffuse  between  cells  and  penetrate  cell  membranes 


2.  OPTOACOUSTIC  PHENOMENA  IN  HETEROGENEOUS  MEDIUM 

Leu  us  for  simplicity  reason  consider  a  turbid  medium  as  a  clear  medium  with  immersed  absorbing  particles  (Rq  -  radius  of 
the  particle).  Let  concentration  of  the  particles,  tu  be  sufficiently  high,  so  that  the  quantity  of  the  particles  in  the 

“elementary”  volume  of  a  x  a  x  a  (see  Fig.l)  is  substantially  greater  then  unity;  »  i  . 


a«a«a 


Figure  1.  Schematics  of  heterogeneous  medium  with  distributed  nanoparticles. 


On  the  other  hand,  the  concentration  should  be  limited,  so  that  the  optical  absorption  over  the  '‘elementary”  volume  can  be 
taken  as  homogeneous:  ana«\  (G  =  ^7tRQ  -  cross-section  of  optical  absorption  of  the  particle).  Particle  absorbs  laser 
radiation  and  it  temperature  rise  up.  Surrounding  liquid  is  heated  due  to  heat  diffusion  from  the  particle.  Pulsed  laser  heating 
of  absorbing  particle  immersed  in  transparent  liquid  can  be  treated  in  several  stages. 

If  the  temperature  of  the  particle  does  not  exceed  the  boiling  point  of  the  liquid,  (lower  limit  of  the  laser  fluence),  the 
phase  transition  cannot  occur. 


Figure  2.  The  system  of  coordinates  for  calculation  of  a  .small  volume  irradiation. 

The  temperature  of  the  liquid  around  the  particle  after  irradiation  by  the  laser  pulse  can  be  taken  in  the  following  form  (the 
temperature  fields  of  the  particles  don’t  overlap  due  to  short  laser  pulse  duration  and  low  concentration  of  the  particles): 

T  =  Tq  +{Tp  -ro)^exp(  -  ~  ^  );  r  >  «o  .  (1) 

where  X  is  diffusivity  of  the  liquid,  is  las^r  pulse  duration.  7„  <  7'^,,,  is  the  temperature  of  the  particle  (we  consider  the 
temperature  distribution  in  the  particle  homogeneous). 

According  to  (1),  the  temperature  of  the  particles  can  be  found  us 
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where  p  is  the  density  of  the  medium  (subscript  /  stands  for  “liquid”,  and  subscript  p  stands  for  “particle”),  is  the  specific 
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heat,  x  =  J X'^T  /  Rq  is  the  ratio  of  heat  diffusion  length  and  particle  radius,  C  = - ^  -  characterizes  the  ratio  of  heat 
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capacity  of  the  particle  and  surrounding  liquid. 

Evaporation  starts  when  the  temperature  of  the  particle  reaches  the  boiling  point.  This  threshold  of  evaporation  can  be 
estimated  as: 


^th  - 


^ley-Jo 


(pCp)[  + 
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The  threshold  has  a  minimum  at  certain  particle  radius,  -tjuin  ~  This  gives  an  optimal  radius  of  the  particle  in  order 

to  minimize  the  evaporation  threshold: 
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(4) 


The  optimal  radius  of  the  particle  depends  on  the  laser  pulse  duration,  specific  heat  of  the  particle  and  specific  heat  of  the 
surrounding  liquid.  For  carbon  or  metal  particle  in  water  and  10  ns  laser  pulse  optimal  radius  is  of  the  order  of  1 10  nm  to  130 
nm.  As  thermal  diffusion  length  in  the  particle  is  of  the  order  of  1  |Lun,  the  model  of  homogeneously  heated  particle  is  valid. 
In  these  conditions  the  threshold  of  evaporation  may  be  presented  as: 


^rnin  =^^^^(pCp)i44]^i4^ +  2^)  ~  14  mJ/cml 


(5) 


In  an  optimum  situation,  heat  containing  in  the  surrounding  liquid  4-5  times  exceeds  that  of  the  particle.  If  the  laser  fluence 
exceeds  the  threshold  (5),  evaporation  at  the  particle  takes  place.  If  laser  fluence  does  not  exceed  the  threshold  (5),  then  the 
evaporation  of  the  liquid  does  not  occur  and  excitation  of  ultrasound  wave  takes  place  due  to  thermal  expansion  of  the 
particle  and  surrounding  layer  of  liquid.  As  the  local  temperature  rise  can  be  several  tens  of  degrees  and  higher,  the  so-called 
thermal  nonlinearity  phenomema  is  manifested.  The  thermal  nonlinearity  leads  to  the  rise  of  the  efficiency  of  optoacoustic 
generation  3-4  times  in  a  case  of  water  (see  below). 

If  the  laser  fluence  exceeds  threshold  (5)  the  temperature  of  the  particle  and  surrounding  liquids  comes  to  boiling  point  at 
some  moment  in  a  course  of  laser  pulse.  Then  thin  vapor  layer  is  generated  around  the  particle.  Due  to  lower  heat 
conductivity  of  vapor  (relatively  to  that  of  liquid),  the  temperature  of  the  particle  continues  to  increase.  After  the  laser  pulse 
the  temperature  of  the  particle  relax  to  the  initial  value  due  to  heat  flux  from  the  particle  through  the  vapor  layer  into  the 
liquid.  The  total  mass  of  the  vapor  is  determined  by  the  heat,  transmitted  from  the  particle  to  the  liquid.  In  a  course  of 
relaxation  at  some  moment  of  time  the  temperature  of  the  particle  becomes  lower  then  the  boiling  point  and  condensation  of 
the  vapor  starts.  Total  mass  of  vapor  produced  around  single  particle  does  not  exceed  the  following  value: 

m^=crOo/A  (6) 

where  X  is  the  heat  of  evaporation.  We  neglect  the  heat  contain  of  the  particle. 

The  pressure  wave,  generated  by  heated  elementary  volume  can  be  taken  in  the  form: 
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where  5V  is  a  variation  of  the  elementary  volume  due  thermal  expansion  and  vapor  production. 

In  case  of  homogeneous  bulk  absorption  of  light  the  expansion  of  an  elementary  volume,  SV ,  can  be  expressed  as  in  [7]: 
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where  is  the  volume  expansion  coefficient.  Thermal  expansion  coefficient  depends  on  the  temperature  and  for  water  this 
dependence  is  quite  strong  (so-called  “thermal  nonlinearity”,  see  Fig. 3).  Therefore,  before  the  temperature  of  the  particle  can 
reach  the  boiling  point,  the  efficiency  of  the  optoacoustic  generation  increases  up  to  2-3  times. 


Figure  3.  Thermal  expansion  coefficient  of  water  as  a  function  of  temperature. 
The  thermal  sources,  producing  the  acoustic  wave  can,  taken  as  follows: 
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Expression  (9)  shows,  that  for  low  temperature  rise 


^  dPjdT  (5r| « 1  the  first  term  in  the  brackets  dominates  and  the 
distribution  of  heat  source  is  proportional  to  the  temperature  rise  .  In  an  opposite  case  of  relatively  high  temperature  rise 
dP/dTSr  »1  the  second  term  is  much  greater  then  the  first  and  the  distribution  of  heat  source  is  proportional  to 

St  .  Thus  for  short  laser  pulse  the  spatial  distribution  of  heat  sources  will  be  a  sum  of  two  exponents: 
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For  low  temperature  rise  (low  fiuence)  the  leading  edge  of  the  optoacoustic  signal  (that  resembles  the  spatial  distribution  of 
heat  sources)  will  have  exponential  shape  with  the  exponent  ;  at  higher  fiuence  the  exponential  slope  will  be  2  times 

sharper:  Ifi^CQt .  In  a  case  of  heterogeneous  medium  with  absorbing  centers,  the  variation  of  the  elementary  volume  due  to 
vapor  production  can  be  estimated  as  (we  consider  pressure  relaxed): 


^2  _  /^a^O 
Pv^ 

where  is  the  vapor  density.  Heating  of  the  vapor  will  give  further  rise  to  the  elementary  volume: 


(11) 


where  yis  the  adiabatic  exponent,  Pq  is  ambient  pressure. 


r-1 

PO  Y 


(12) 


Let’s  compare  the  efficiency  of  various  mechanisms  of  the  optoacoustic  phenomena.  Volume  deformation  per  unit  specific 
heat  release  for  various  mechanisms  are  presented  in  the  Table  below. 


Thermal  expansion 
of  liquid  water 

Vapor 

Production 

Thermal  expansion 
of  water  vapor 

crri^  U 

8.6*10"^ 

3.2*10"^ 

1.7 

Thus,  light  absorption  by  submicron- size  particles  can  enhance  the  efficiency  of  optoacoustic  generation  by  several  orders  of 
magnitude,  but  this  process  has  a  threshold  dependent  on  laser  pulse  duration,  size  of  the  particle  and  the  light  absorption 
cross-section. 


3.  MATERIALS  AND  EXPERIMENTAL  SET-UP 

Studies  of  the  optoacoustic  phenomena  in  heterogeneous  media  were  performed  with  the  Q-s witched  Nd-YAG  laser  (laser 
pulse  duration  ~\4  ns  at  1/e  level)  operated  at  fundamental  (A  =  L064)jLm)  wavelength  and  the  second  harmonic 
( A  =0.532p.m).  The  laser  pulse  energy  was  varied  by  neutral  density  filters.  Pulse  repetition  rate  did  not  exceed  3  Hz  and 
was  chosen  to  fulfill  the  safety  requirement  not  to  elevate  the  average  temperature  of  the  medium  more  than  1  degree. 
Actually,  most  of  the  measurements  have  been  done  with  single  laser  pulses. 

We  used  optically  heterogeneous  media:  water  solution  of  India  ink  of  various  concentrations,  water  suspension  of  gold 
nanoparticles  (diameter  10  nm,  30  nm  and  100  nm)  with  the  volume  concentration  of  the  particles  of  the  order  of  10'^.  As  a 
reference  homogeneous  medium  water  solution  of  cupric  chloride  of  various  concentrations  was  employed. 

The  enhancement  of  the  optoacoustic  imaging  contrast  in  heterogeneous  medium  was  investigated  with  laser  optoacoustic 
imaging  system  (LOIS-2)  designed  for  breast  cancer  imaging.  The  system  uses  32-element  arc-shaped  array  of  piezoelectric 
transducers  with  ultrasonic  frequency  band  of  0.05-3  MHz,  which  provides  spatial  resolution  of  the  order  of  0.5  mm. 


4.  EXPERIMENTAL  RESULTS 

The  leading  edges  of  optoacoustic  signals  generated  in  1%  solution  of  India  ink  in  water  for  various  energy  fluences  are 
depicted  in  Fig.  4.  The  light  absorption  coefficient  of  this  solution  was  =  18  cm~^ . 


Figure  4.  Profiles  of  leading  edges  of  the  optoacoustic  signals  in  aqueous  solutions  of  1%  India  ink  generated  by  laser  pulses 
with  three  different  fluences. 
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For  low  fluence  (Oq  =14.7  mJ  lent  ,  curve  I)  the  leading  edge  has  exponential  shape  with  the  exponent,  consistent  with 

the  light  absorption  coefficient  (see  (10)).  With  an  increase  of  the  incident  fluence  (  3>q  =  88.7  mJ  Icm^y  curve  3),  the  slope 
becomes  steeper  and  corresponds  to  2  times  sharper  exponent  (see  (10),  the  second  term).  In  the  intermediate  case 

($0  =  32.4  mJ  /  cm  ,  curve  2)  beneath  the  surface  the  slope  is  two  times  higher,  then  in  the  depth  of  the  irradiated  volume 
(see  Fig.4).  This  is  in  consistence  with  the  “thermal  nonlinearity”  model  (see  (9)-(10)). 


Let’s  discuss  the  intermediate  case  Oq  =  32.4  mJ  I  cm  in  more  details.  If  we  assume  homogeneous  absorption  of  light  in 
the  aqueous  solution,  the  maximum  temperature-rise  will  be  AT"  =/i^Oo/pc^  s;0.14  On  the  other  hand,  the  influence 

of  thermal  nonlinearity  is  so  significant,  that  temperature  elevation  can’t  be  expected  such  low.  The  local  temperature-rise 
can  be  estimated  from  the  ratio  of  the  amplitudes  of  exponents  in  (10).  Fit  of  the  leading  edge  of  the  optoacoustic  signal  with 

1 

the  formula  (10)  yields  Pq  - ST  -0.51,  which  makes  it  possible  to  estimate  the  local  temperature  elevation  as 

dT 


5T  ^\0K  (for  water  at  room  temperature  dp/dT  -  0.05  K  ^ ).  Thus,  the  absorbing  carbon  particles  are  exposed  to 
at  least  70  times  higher  temperature-rise  then  the  average  temperature  elevation  in  the  irradiated  volume  of  the  solution. 


Flux  (mJ/errf) 


Figure  5.  Amplitude  of  optoacoustic  signal  vs.  energy  fluence  for  water  suspension  of  carbon  particles. 


The  dependence  of  the  amplitude  of  optoacoustic  signal  vs.  laser  energy  fluence  for  water  suspension  of  carbon  particles  is 
presented  in  Fig.5.  The  diameter  of  the  particles  was  approximately  0.5  pm  and  was  close  to  optimal  value  (see  (4)),  so  it 
was  reasonable  to  expect  manifestation  of  evaporation  at  absorbing  centers  for  the  laser  fluence  of  approximately 
2 

(15  -  20)  mJ  /  cm  and  significant  thermal  nonlinearity  below  this  limit. 


The  presented  data  is  consistent  with  the  theoretical  model  discussed  above.  For  fluences  Oq  <  8  mJ I cm^  the  amplitude  of 
the  optoacoustic  signal  generated  in  heterogeneous  medium  (points  in  Fig.5)  is  proportional  to  the  fluence  and  coincides  well 
with  the  theoretical  value  for  homogeneous  medium  (solid  line).  This  means  that  the  temperature-rise  at  the  absorbing 
particle  is  relatively  low  (at  least  below  10  K)  and  no  thermal  nonlinearity  is  manifested. 

2  2 

In  the  fluence  range  of  8  mJ  / cm  <  Oq  <15  mJ  / cm  the  amplitude  of  the  optoacoustic  signal  generated  in  heterogeneous 
medium  exceeds  that  in  a  homogeneous  medium  up  to  3  times.  This  corresponds  well  to  the  rise  of  optoacoustic  generation 
efficiency  due  to  thermal  nonlinearity  below  .the  evaporation  threshold.  For  higher  fluence  Oq  >15  a  sharp 

increase  of  the  amplitude  of  excited  signals  takes  place.  It  may  occur  due  to  the  elevation  of  the  temperature  of  particles  to 
the  boiling  point  and  following  evaporation  of  water  layer  around  the  absorbing  particle.  In  accordance  with  the  enhancement 
of  the  efficiency  of  optoacoustic  generation  in  vapor  production  regime  (see  Table),  a  dramatic  increase  of  the  amplitude  of 


optoacoustic  signal  can  be  seen.  The  enhancement  of  the  amplitude  of  optoacoustic  signal  about  2  orders  was  achieved.  The 
threshold  of  evaporation  at  the  particle  is  in  good  agreement  with  the  theoretical  estimations  (see  (5)). 

The  experiments  described  above  were  performed  with  a  single  piezoelectric  transducer.  To  investigate  the  enhancement  of 
the  contrast  of  optoacoustic  imaging  with  the  heterogeneous  absorbing  medium,  a  32-elements  arc-shaped  transducer  array  of 
LOIS-2  imaging  system  was  employed.  The  transducer  was  attached  to  the  tank,  filled  with  milk  (light  attenuation  coefficient 

j^ejf  Two  2.4-mm  diameter  vessels  were  embedded  in  the  tank.  One  of  them  was  filled  with  water  solution  of 

cupric  chloride  (homogeneous  absorber),  the  other  -  with  water  solution  of  India  ink  (heterogeneous  absorber).  The  light 
absorption  coefficient  of  heterogeneous  and  homogeneous  medium  was  the  same. 
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Figure  6.  Optoacoustic  image  of  vessels  with  heterogeneous  (7)  and  homogeneous  (2)  absorption. 
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Figure  7.  Optoacoustic  image  of  vessels  with  gold  particles  10  nm  in  diameter  (7),  30  nm  (2)  and  100  nm  (3). 

2  2 
Oq  =  13.4  mJ  / cm  (left  panel).  Optoacoustic  image  of  similar  nanoparticles  at  Oq  =71  mJ  I (right  panel). 


Optoacoustic  images  of  the  same  vessels  filled  with  the  water  suspension  of  gold  particles  of  various  diameter  at 
X  =  0.53 //m  are  presented  in  Fig.7  (volume  concentration  of  the  particles  was  3*10“^).  For  the  fluence  below  the  evaporation 
threshold  (Fig.7a)  the  100  nm  in  diameter  particles  (vessel  3)  produce  the  brightest  image,  due  to  higher  light  absorption 
coefficient  of  the  solution.  Correspondingly,  the  solutions  of  30  nm  (vessel  2)  and  10  nm  (vessel  1)  particles  produce  the 
images  with  lower  intensity.  For  the  fluence  above  evaporation  threshold  (Fig.7b)  the  brightest  image  was  produced  by  10 
mri  particles,  and  the  intensity  of  the  images  is  reverse  to  that  below  the  threshold. 


5,  DISCUSSION  AND  CONCLUSIONS 

Light  absorbing  nanoparticles  were  studied  as  a  potentially  very  effective  contrast  agent  for  laser  optoacoustic  imaging  in 
tissue.  Enhanced  efficiency  of  acoustic  wave  generation  can  take  place  due  to  superheating  of  absorbing  particles  followed 
by  evaporation  of  the  surrounding  water  layer.  Some  other  phenomena  of  phase  transition  and  gas  production  can  be  utilized. 
Laser  generation  of  ultrasound  in  heterogeneous  medium  can  be  2-3  orders  more  efficient  than  that  in  optically  homogeneous 
medium.  An  optimal  particle  size  yielding  minimal  threshold  of  evaporation  at  the  particle  surface  was  found  in  the  range  of 
~200-nm.  The  threshold  of  evaporation  is  relatively  high  and  exceeds  10  mJ/cm^  for  metal  or  carbon  particles  suspension  in 
water.  However,  application  of  the  plasmon  resonance  absorption  in  the  particle,  the  threshold  can  be  lowered  drastically 
(due  to  higher  crossection  of  light  absorption).  One  can  expect  a  decrease  of  the  threshold  up  to  2-3  orders  of  magnitude  in 
case  of  utilization  of  plasmon  resonance  for  nanoparticle  absorption. 

Nanotechnology  is  a  rapidly  expanding  field  of  science  and  engineering  that  offers  a  unique  potential  for  integrating  contrast 
agents  and  receptor-specific  markers  for  cancer,  promising  to  yield  smart  contrast  agents.  Gold  nanoshells  filled  with  silica 
can  be  produced  to  absorb  strongly  the  desired  wavelength  of  near  infrared  laser  pulses  [14],  These  nanoparticles  may  be 
further  modified  to  yield  exponential  enhancement  of  optoacoustic  signals  upon  pulsed  laser  irradiation.  Effective 
bioconjugates  have  been  found  to  link  various  nanoparticles  with  biological  molecules  [15].  The  problem  of  selective 
targeting  of  cancerous  cells  and  tissues  is  currently  also  under  active  development.  Peptides,  other  ligands  and  antibodies 
specific  to  surface  receptors  in  breast  cancer  cells  or  endothelial  cell  of  tumor  microcapillary  network  are  available  for 
selective  targeting  of  malignant  tumors  [16,17].  Thus,  all  components  are  available  for  successful  development  of  the 
optoacoustic  imaging  technique  with  supercontrast.  A  receptor-specific  contrast  agent  in  combination  with  optoacoustic 
imaging  may  dramatically  enhance  the  level  of  sensitivity  of  detection  of  small  breast  tumors,  reduce  the  rate  of  false 
positives  and  negatives,  provide  diagnostic  imaging  information,  and  be  useful  in  monitoring  for  early  stages  of  breast  cancer 
and  cancer  recurrence. 

We  conclude  that  absorbing  nanoparticles  conjugated  with  targeting  agents  specific  to  cancer  may  be  applied  as  a  particulate 
contrast  agent  for  early  detection  of  cancer  with  laser  optoacoustic  imaging  system. 


6.  ACKNOWLEDGEMENTS 

This  work  was  supported  by  the  National  Cancer  Institute  (grant  #R29-CA80221)  and  DOD  Breast  Cancer  Research 
Program,  US  Army  (grant  #DAMD  17-99- 1-9404). 


7.  REFERENCES 

1.  A.A.  Oraevsky,  S.L.  Jacques,  R.O.  Esenaliev,  F.K.  Tittel:  Time-Resolved  Optoacoustic  Imaging  in  Layered 
Biological  Tissues",  In:  "Advances  in  Optical  Imaging  and  Photon  Migration",  vol.  21,  ed.  by  R.R.  Alfano, 
Academic  Press  (1994)  pp.  161-165. 

2.  A.A.  Oraevsky,  S.L.  Jacques,  R.O.  Esenaliev,  F.K.Tittel:  Laser  based  optoacoustic  imaging  in  biological  tissues, 
Proc.  SPIE  1994;  2134A:  122-128. 

3.  R.A.  Kruger,  P.  Liu:  Photoacoustic  ultrasound:  Pulse  production  and  detection  in  0.5%  Liposyn,  Medical  Physics, 
1994;  21(7):  1179-1184, 

4.  R.A.  Kruger:  Photoacoustic  ultrasound,  Med,  Phys.  1994;  21(1):  127-131 

5.  A.A.  Oraevsky,  S.L.  Jacques,  F.K.  Tittel:  Measurement  of  tissue  optical  properties  by  time-resolved  detection  of 
laser-induced  transient  stress,  Applied  Optics,  1997;  36(1):  402-415. 


6.  A.A.  Oraevsky:  Laser  optoacoustic  imaging  for  cancer  diagnosis,  LEOS  Newsletter  1996;  10(6):  17-20. 

7 .  A.A.  Oraevsky,  V.G. Andreev,  A.A.Karabutov,  R.O.Esenaliev,  “Two-dimensional  optoacoustic  tomography: 
transducer  array  and  image  reconstruction  algorithm”,  SPIE  Proceed,  3601,  256-267  (1999). 

8.  R.O.  Esenaliev,  A.A.  Karabutov,  A.A.  Oraevsky:  Sensitivity  of  laser  opto-acoustic  imaging  in  detection  of  small 
deeply  embedded  tumors,  lEEEJ.  ST  Quant.  Electr.  1999;  5(4):981-988. 

9.  G.J.  Diebold,  M.L  Khan,  S.M.  Park:  “Photoacoustic  signatures  of  particulate  matter:  optical  production  of  acountis 
monopole  radiation”,  Science  1990;  250:  101-104. 

10.  F.V.  Bunkin,  M.L  Tribelsky,  Sov.Phys.Uspekhi,  130,  193  (1980). 

11.  M.W.  Sigrist,  Laser  generation  of  acoustic  pulses,  J.  Appl.  Phys.,  60,  R83  (1986). 

12.  R.O.  Esenaliev,  A.A.  Karabutov,  N.B.  Podymova,  V.S.  Letokhov,  “Laser  ablation  of  aqueous  solutions  with 
spatially  homogeneous  and  heterogeneous  absorption”,  Appl.  Phys.  B,  59,  73-81  (1994). 

13.  A.C,  Beveridge,  T.E.  McGrath,  G.J.  Diebold,  A.A.  Karabutov,  “Photoacoustic  shock  generation  in  carbon 
suspension”,  Appl.  Phys.  Lett.,  75(26),  4204-4206  (1999). 

14.  S.J.  Oldenburg,  R.D.  Averitt,  S.L.  Wescott,  N.J.  Halas:  Nanoengineering  of  optical  resonances,  Chem.  Phys.  Lett. 
1998;  288:  243-247. 

15.  H.  Pinto- Alphandary,  A.  Adremont,  P.  Couvreur:  “Targeted  delivery  of  antibiotics  using  liposomes  and 
nanoparticles:  research  and  applications.  Int  J.  Antimicrob  Agents  2000;  13(3):155-168. 

16.  H.M.  Ellerby,  R.  Pasqualini  et  al.  Anti-cancer  activity  of  targeted  pro-apoptotic  peptides,  Nature  Med.  1999;  5(9): 
1032-1038. 

17.  L.M.  Weiner,  G.P.  Adams:  New  approaches  to  antibody  therapy,  Oncogene  2000;  19(53):  6144-6151. 


Optoacoustic  images  of  early  cancer  in  forward  and  backward  modes 

A.A.  Karabutov*,  V.G.  Ancireev\  B.  Bell',  R.D.  Fleming\  Z.  Gatalica^  M.  Motamedi' ^ 

E.V.  Savateeva',  H.  Singh^  S.V.  Solomatin',  S.L.  Thomsen^  P.M.  Henric’hs^  A.A.  Oraevs’ky' 

'Center  for  Biomedical  Engineering,  ^Department  of  Pathology,  ^Department  of  Radiology, 

^  "^Department  of  Surgery,  University  of  Texas  Medical  Branch  at  Galveston,  Texas 

^LaserSonix  Technologies  Inc.,  Houston,  Texas,  ^MD  Anderson  Cancer  Center,  Houston,  Texas 


ABSTRACT 

Optoacoustic  tomography  combines  advantages  of  pronounced  optical  contrast  between  different  tissues  and  high  resolution 
of  wide-band  ultrasound  imaging.  Laser  pulses  may  be  effectively  used  to  produce  acoustic  sources  in  tissue  with  enhanced 
optical  absorption.  Ultrasonic  waves  can  propagate  in  biological  tissue  with  minimal  distortion  and  deliver  diagnostic 
information  to  the  surface  of  tissue,  where  they  may  be  detected  with  wide-band  ultrasonic  transducers.  Current  status  of  the 
optoacoustic  tomography  applied  in  early  detection  of  cancerous  lesions  in  the  breast  (utilizing  forward  mode)  and  in  oral 
cavity  (utilizing  backward  mode)  is  reviewed. 


Keywords.  Optoacoustic  imaging,  in  vivo  diagnostics,  breast  cancer,  ultrawide- band  acoustic  transducer,  laser,  ultrasound 


1.  INTRODUCTION 

Laser  optoacoustic  imaging  system  (LOIS)  was  proposed  for  cancer  detection  in  order  to  combine  advantages  of  optical 
contrast  and  sensitive,  high-resolution  ultrasonic  detection  in  one  tomography  system  [1-4].  One  of  the  major  applications  of 
LOIS  is  the  detection  of  early  breast  cancer  [5-7].  Other  cancerous  lesions  (such  as  in  prostate  and  digestive  system)  having 
optical  properties  different  from  normal  tissues  also  can  be  detected  with  optoacoustic  tomography  (OAT)  [8-10].  OAT 
utilizes  analysis  of  profiles  of  acoustic  (ultrasonic)  signals  induced  by  laser  pulses  in  tissue  with  tumors  [11-13].  Application 
of  the  ultrawide-band  ultrasonic  detection  instead  of  detection  of  photons  helps  to  overcome  problems  associated  with  strong 
optical  scattering  in  biological  tissues  and  improve  depth  of  monitoring,  sensitivity  and  spatial  resolution  [14-16].  The 
profiles  of  acoustic  waves  generated  under  irradiation  conditions  of  temporal  pressure  confinement  in  the  volume  of  tissue 
resemble  the  profiles  of  absorbed  laser  energy  in  the  tissue,  and  can  be  used  to  replicate  the  tissue  structure.  Tumors  with 
dimensions  of  1.5-mm  to  15-mm  irradiated  with  laser  pulses  represent  themselves  as  sources  of  acoustic  waves  with 
ultrasonic  frequencies  of  -1  MHz  to  -100  kHz.  Such  ultrasonic  waves  can  propagate  in  biological  tissues  with  insignificant 
attenuation  [17,18].  Therefore,  tumors  in  the  breast  can  be  detected  at  distances  of  up  to  7-8  cm.  However,  spherical 
propagation  of  acoustic  waves  changes  initial  profile  of  absorbed  laser  energy  (especially  in  the  ultrasonic  frequency  range  < 
1  MHz),  so  that  the  detected  signals  must  be  processed  by  integrating  them  over  the  entire  time-course  of  detection  in  order 
to  be  used  in  the  optoacoustic  image  reconstruction  [5-7].  Smaller  cancerous  lesions  or  other  tissue  structures  produce  wider 
range  of  ultrasonic  frequencies  up  to  100  MHz  [17].  Due  to  stronger  attenuation  of  high-frequency  ultrasound,  optoacoustic 
profiles  arising  from  microscopic  lesions  and  other  profiles  with  fine  structure  or  can  be  picked  up  from  shallower  depths  of 
<  1-cm.  The  detection  of  ultrasonic  waves  can  be  achieved  with  ultrawide-band  piezoelectric  transducers  [18-20]  or  optical 
interferometers  [21-23].  The  advantage  of  the  former  is  superior  sensitivity  and  low  noise  per  unit  bandwidth.  The 
advantage  of  the  latter  is  noncontact  detection  and  rapid  examination  of  large  area. 

There  are  two  modes  of  detection  for  the  laser-induced  pressure  transients  (LIFT),  termed  as  “forward”  and  “backward” 
modes.  In  both  cases  the  irradiation  of  the  surface  of  tissue  under  evaluation  takes  place  through  an  optically  and  acoustically 
transparent  buffer  (such  as  quartz,  water  or  plastic).  This  “optoacoustic  buffer”  permits  good  acoustic  coupling  to  the  tissue 
and  incorporation  of  light  delivery  system  and  detection  system.  Usually  as  wide  laser  beam  is  employed  for  tissue 
illumination.  The  thermal  pressure  arising  inside  heated  volume,  due  to  rapid  deposition  of  laser  energy,  yields  ultrasonic 
transients  launched  in  two  opposite  directions,  one  into  absorbing  and  the  other  one  into  transparent  buffer.  The  temporal 
shape  of  the  LIPTs  depends  on  the  ratio  of  acoustic  impedances  of  tissue  and  optoacoustic  buffer  [16-18]  and  may  be 
different  for  LIPTs  propagating  forward  into  tissue  and  backward  into  optoacoustic  buffer.  In  the  forward  mode,  the 
detection  is  being  made  of  LIFT  launched  into  the  tissue  forward  along  the  direction  of  the  laser  beam.  The  ultrasonic 
detection  takes  place  at  the  rear  surface  of  the  irradiated  tissue. 


The  temporal  shapes  of  LEUT  propagated  into  the  absorbing  medium  are  well  established  and  quantitatively  investigated 
both  theoretically  and  experimentally  (see,  for  example  [18,24-26]).  Their  use  is  perspective  for  numerous  applications,  but 
in  a  case  of  the  forward  mode  detection  one  has  to  have  a  two-sides  access  to  the  medium  under  investigation  -  both 
irradiated  and  rare  surfaces  of  the  medium  have  to  be  accessible.  The  backward  mode  of  LEUT  detection  eliminates  this 
problem.  In  this  case  the  detection  of  ultrasonic  transient  takes  place  at  the  rear  surface  of  transparent  medium,  coupled  to 
investigated  one  (the  irradiation  of  the  absorbing  medium  is  realized  through  the  side  of  transparent  medium).  As  the 
transparent  medium  can  be  taken  smooth  enough,  an  additional  ultrasonic  scattering  and  absorption  in  a  course  of  ultrasonic 
wave  propagation  can  be  eliminated. 

The  two  detection  modes  allow  certain  flexibility  for  in  vivo  and  in  vitro  measurements  and  enables  one  to  investigate  both, 
the  absorbed  optical  energy  distribution  [22,23]  and  ultrasonic  back-scattering  [18,26].  This  paper  reviews  optoacoustic 
tomography  of  the  breast  cancer  in  the  forward  mode  and  imaging  of  early  cancer  in  orla  mucosa  using  the  backward  mode. 

The  initial  years  of  our  research  in  the  area  of  optoacoustic  imaging  of  the  breast  yielded  the  following  results:  (1)  first 
clinical  prototype  system  for  the  breast  imaging  was  developed,  fabricated  and  extensively  tested  in  phantoms  [10,12],  (2) 
detection  and  localization  of  breast  tumors  in  surgical  mastectomy  specimens  were  demonstrated  [15],  (3)  parameters  of 
LOIS  were  compared  with  X-ray  mammography,  MRI  and  ultrasound  imaging  [15].  The  main  merit  of  optoacoustic  imaging 
was  found  in  the  possibility  to  utilize  high  tissue  contrast  based  primarily  on  optical  absorption  and  simultaneously  preserve 
high  resolution  of  ultrasonic  imaging.  The  image  contrast  is  defined  as  the  difference  in  optoacoustic  image  brightness  in 
tumors  and  in  adjacent  normal  tissue  divided  to  the  brightness  of  normal  tissue  background.  The  value  of  optoacoustic 
contrast  was  found  in  the  range  of  1  to  5,  which  substantially  exceeds  any  other  endogenous  tissue  contrast  currently  utilized 
in  medical  imaging:  ultrasonography,  MRI  and  X-ray  radiography.  Based  on  preliminary  comparative  analysis  of  gross 
tumor  cross-sections,  immunohistology  and  optoacoustic  images  we  hypothesize  that  the  optoacoustic  contrast  results 
primarily  from  increased  optical  absorption  in  the  dense  micro  vascularity  of  the  tumors  [7,12].  However,  optical  scattering 
in  tumors,  especially  those  treated  with  chemo-  and  radiotherapy,  may  also  significantly  contribute  to  the  optoacoustic 
contrast  [27].  Thermoelastic  properties  of  the  breast  tumors  may  also  add  to  the  optoacoustic  contrast,  however,  were  not  yet 
quantified  [28].  Comparison  of  optoacoustic  images  with  those  of  X-ray  radiography,  MRI  and  ultrasound  imaging 
performed  in  vitro  revealed  good  correlation  in  tumor  size  and  location,  and  the  spatial  resolution  equal  to  that  of  ultrasound 
(better  than  1  mm)  was  achieved  [15]. 

In  continuous  effort  of  developing  an  imaging  modality  with  better  spatial  resolution  applicable  for  sensitive  detection  of 
variety  of  cancerous  lesions,  our  group  designed  the  optoacoustic  front  surface  transducers  (OAFST)  operating  in  the 
backward  mode,  i.e.  detecting  acoustic  profiles  at  the  site  of  laser  irradiation  [29].  Feasibility  study  in  characterization  of 
tissue  structure  performed  with  of  OAFST  confirmed  its  utility  for  in  vivo  imaging  of  various  skin  lesions  and  cancer  staging 
in  oral  cavity  [30,31].  An  ultrawide  band  of  ultrasonic  detection  realized  in  OAFST  yields  in-depth  resolution  of  about  15- 
jim.  In  order  to  achieve  lateral  resolution  of  OAFST  comparable  with  in-depth  resolution,  a  focused  optoacoustic  front 
surface  transducer  was  developed  [28].  The  new  transducer  design  is  referred  to  as  Confocal  Opto- Acoustic  Transducer 
(COAT).  The  performance  of  COAT  was  initially  tested  in  optically  homogeneous  and  heterogeneous  phantoms  [29,31]. 
Experiments  in  vivo  also  were  performed  demonstrating  capability  of  the  developed  transducer  in  imaging  and  distinguishing 
early  stages  of  squamous  cell  carcinoma  in  oral  mucous  of  golden  hamsters  [23-25]. 

2.  OPTOACOUSTIC  TOMOGRAPHY  OF  DEEP  TUMORS  IN  THE  BREAST 

2.1.  Laser  Optoacoustic  Imaging  System 

Schematic  diagram  of  LOIS  is  shown  in  Fig.  la.  An  Nd:YAG  laser  (Big  Sky  Lasers,  MA)  operating  at  1064  nm  was  used  as 
a  source  of  pulses  of  10-ns  duration  with  repetition  rate  of  20  Hz.  A  quartz  optical  fiber  was  employed  for  the  laser  pulse 
delivery  to  the  tissue  surface.  A  telescope  was  used  to  expand  laser  beam  from  the  fiber  and  deliver  a  parallel  beam  with 
diameter  of  1-cm  to  the  surface  of  the  skin.  A  32-element  arc-shaped  array  of  PVDF  acoustic  transducers  was  employed  in 
LOIS  for  detection  of  LIPTs  within  ultrasonic  frequency  band  of  20  kHz  to  2  MHz  (see  Fig.l).  An  aperture  angle  of  120 
degrees  provides  lateral  resolution  of  about  1.2-mm,  comparable  with  in-depth  resolution  of  0.8 -mm,  determined  by 
frequency  band  of  the  transducer.  On  the  other  hand,  this  geometry  provides  sufficient  resolution  within  the  field  of  view  of 
50-mm  X  80-mm  around  the  focal  zone.  The  Sensitivity  of  the  piezoelectric  transducers  was  S  =  (10±2)  juV/Pa.  The  noise 
level  of  the  transducer  with  the  preamplifier  corresponded  to  the  pressure  value  2.5-3  Pa.  The  dynamic  range  of  the  ultrasonic 
detection  system  was  estimates  as  73  dB.  These  specifications  make  it  possible  to  obtain  high  quality  optoacoustic  image  at 
the  depth  of  up  to  80-mm  using  laser  fluence  of  20  mJ/cm^  incident  on  the  irradiated  surface. 


Figure  1.  Schematic  diagram  of  clinical  LOIS  (a)  and  a  photograph  of  the  32-element  arc  transducer  array  (b). 

Various  band-pass  filters  were  applied  for  optoacoustic  signal  treatment.  They  were  used  for  elimination  of  exponential 
background,  arisen  due  to  regular  attenuation  and  absorption  of  laser  radiation  in  tissue.  The  location  of  irradiated  surface 
area  relatively  to  the  transducer  was  determined  from  the  optoacoustic  signals  and  marked  with  bright  line. 

A  radial  back-projection  algorithm,  developed  for  the  optoacoustic  imaging  and  described  in  [10,  32-34]  was  employed.  For 
the  reconstruction  of  two-dimensional  optoacoustic  images  we  used  signal  integrals  projected  back  onto  the  two-dimensional 
grid  taking  into  account  the  directivity  pattern  (angle  of  acceptance)  of  each  transducer  in  the  array.  The  image  represents 
distribution  of  a  product  of  thermo- acoustic  efficiency,  optical  absorption  and  absorbed  laser  energy.  The  back-projected 
images  acquired  with  only  limited  number  of  detectors  display  artifacts  associated  with  and  image  reconstruction  using 
incomplete  set  of  data.  The  contrast  of  these  images  could  be  improved  through  the  image  filtration  procedure  applied  to  the 
entire  two-dimensional  image  [15].  However,  unprocessed  optoacoustic  images  contain  quantitative  information  of 
distribution  of  absorbed  optical  energy  in  tissue.  This  quantitative  information  may  be  used  for  functional  imaging  and 
estimate  of  tissue  optical  properties. 

Data  acquisition  in  LOIS  is  realized  with  high  efficient  32-channel  12-bit  ADC.  The  laser  pulse  repetition  rate  was  20  Hz 
and  averaging  over  16  pulses  are  employed.  Data  processing  and  back- projected  image  formation  takes  0.8  -3  seconds 
(depending  on  the  quantity  of  pixels  in  image).  Image  filtration  takes  0.2  seconds.  So,  total  time  equals  2-4  sec.  Current 
design  of  the  laser  optoacoustic  imaging  system  (LOIS)  allows  operator  to  scan  the  acoustic  detector  array  manually  through 
the  breast  and  to  obtain  images  of  chosen  parts  of  the  breast. 

2.2.  Initial  Tests  of  Clinical  LOIS 

The  is  a  strong  evidence  that  malignant  breast  tumors  have  significantly  enhanced  content  of  deoxygenated  blood  in  the 
tumor  angiogenesis  [35,36].  We  hope  that  based  on  imaging  hemoglobin  vs  deoxyhemo  spectroscopic  optoacoustic 
Therefore,  the  Laser  Optoacoustic  Imaging  system  was  tested  in  phantoms  resembling  breast  with  blood  vessels.  As  shown  in 
the  diagram  on  Figure  2,  the  artificial  blood  vessels  were  located  at  different  depths  and  were  placed  either  parallel  or 
perpendicular  to  the  imaging  plane.  Position  of  the  transducer  array  was  at  the  bottom  of  the  image  (not  shown).  Surface 
position  is  visible  as  a  bright  area  on  the  images.  A  single  fiber  was  used  for  irradiation  at  only  one  location  on  the  surface  of 
the  phantom.  A  photograph  of  the  phantom  is  presented  in  Fig.  2  (right).  The  phantom  was  imaged  twice,  one  time  with 
optical  fiber  positioned  at  one  point  of  its  surface  and  the  second  time  at  the  opposite  position  (the  phantom  was  rotated  1 80 
degrees). 

Two-dimensional  images  of  blood  vessels  of  different  diameter  and  filled  with  blood  having  different  level  of  oxygenation  is 
shown  in  Figure  3.  Whole  milk  was  diluted  by  water  to  obtain  the  effective  optical  attenuation  coefficient  1.4  cm'^  in  order  to 
mimic  the  upper  limit  of  optical  scattering  in  the  breast  tissue  at  the  wavelength  of  1064  nm.  Optical  absorption  coefficients 
of  blood  in  vessels  varied  from  0.8  to  4.0  per  cm  depending  on  oxygen  content  The  first  image  was  reconstructed  using 
filtered  optoacoustic  signal  in  order  to  remove  the  exponential  trend  associated  with  effective  optical  attenuation  (Fig.  3  left). 
The  second  image  was  not  filtered,  so  that  exponential  profile  of  light  distribution  in  the  phantom  could  be  observed  (Fig.  3, 
right). 


1.6  mm  dia 
2.5  cm  depth 
17%  O2 


Figure  2.  Left  panel:  Collagen  gel  phantom  with  artificial  blood  vessels  filled  with  rabbit  blood  having  various 
levels  of  oxygenation.  Positions  of  blood  vessels  and  percent  of  oxyhemoglobin  are  shown  in  windows.  Optical 
properties  of  this  phantom  were  similar  to  the  optical  properties  of  the  breast  at  the  wavelength  of  1064  nm.  Right 
panel:  Photograph  of  the  same  phantom  with  artificial  blood  vessels.  Transducer  array  and  optical  fiber  are  also 
presented  on  the  photograph.. 
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Figure  3.  Optoacoustic  image  of  blood  vessels  in  milky  gelatin  phantom.  Contrast  of  blood  vessels  relative  to 
background  is  shown  on  Contrast  Bar.  Depth  from  the  surface  is  depicted  on  the  left  vertical  axis  of  each  image. 

The  image  of  blood  vessel  cross-sections  differs  from  the  correct  circular  shape.  This  fact  can  be  explained  by  limited  angle 
aperture  of  the  transducer.  This  limited  lateral  resolution  in  the  far  zone  of  the  images  presented  in  Fig.  3.  An  increase  of  the 
acceptance  angle  for  the  blood  vessel  located  farther  from  the  irradiation  site  and  closer  to  the  detector  array  resulted  in  a 
better  reproduction  of  its  circular  cross-section  in  the  image,  as  clearly  depicted  in  Fig.  3  (right  panel)  showing  an  image  of  a 
blood  vessel  closer  to  the  transducer.  The  angle  of  acceptance  in  this  case  is  about  180°  and  as  a  consequence  the  shape  of 
image  highly  correlates  with  real  shape  of  the  circle.  Note  that  our  image  reconstruction  algorithm  allows  accurate 
reproduction  of  spherical  objects  and  cylindrical  objects  (see  a  blood  vessel  in  the  center  going  in  horizontal  direction). 
Thus,  we  may  predict  that  objects  with  random  shapes  will  be  correctly  reproduced  on  LOIS  images. 
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Figure  4.  Optoacoustic  image  of  human  arm  in  the  area  of  byceps  (left  panel)  and  m-depth  cross-section  of  this  image  resembling  the 
profile  of  absorbed  optical  energy  (right  panel).  The  effective  optical  attenuation  coefficient,  =1.28  cm'b 


All  blood  vessels  were  clearly  resolved  and  depicted  as  separate  objects.  The  relative  position  of  the  blood  vessels  as 
depicted  on  the  optoacoustic  images  accurately  resembles  their  position  as  shown  on  photograph  Fig.  2  (right  panel)  and  their 
depths  are  in  good  agreement  with  their  real  location.  This  imaging  experiment  demonstrated  sensitivity  of  LOIS  sufficient 
for  detection  of  blood  vessels  at  the  depth  of  over  75  mm  in  an  optical  scattering  and  absorbing  phantom.  The  brightness  of 
the  blood  vessel  images  was  proprtional  to  the  amount  of  absorbed  optical  energy,  which  in  turn  was  proportional  to  the  level 
of  blood  oxygenation  in  each  vessel. 

2.3.  Quantitative  Optoacoustic  Imaging 


Optoacoustic  tomography  visualizes  the  spatial  distribution  of  absorbed  optical  energy  in  tissue.  However,  distortion  of 
optoacoustic  signals  upon  propagation  through  thick  layers  of  tissue  and  signal  processing  may  in  principle  modify  the 
original  profile.  We  tested  possibility  to  obtain  axial  in  depth  profiles  of  absorbed  optical  energy  from  optoacoustic  images 
from  human  muscle  tissue.  A  human  arm  was  imaged  in  the  area  of  byceps  with  LOIS  designed  for  breast  cancer  imaging 
(see  Fig. 4). 


Only  high  frequency  filtering  was  applied  to  the  optoacoustic  image.  This  experiments  demonstrated  that  quantitative 
information  can  be  obtained  form  optoacoustic  images.  Neither  signal  correction  for  acoustic  diffraction,  nor  reconstruction 
of  absorbed  optical  energy  distribution  from  optoacoustic  signals,  nor  correction  for  sensitivity  of  individual  piezoelectric 
transducers  changed  quantitative  information.  The  effective  optical  attenuation  coefficient,  |L4ff=1.28  cm'^  determined  form 
exponentil  fit  of  the  optoacoustic  profile,  was  in  a  good  agreement  with  results  of  direct  measurements  of  optical  properties  in 
this  tissue. 

Imaging  Breast  Cancer  in  vivo 

Eligible  patients  were  chosen  from  the  group  of  patients  scheduled  for  radical  surgical  mastectomy  with  breast  cancer 
diagnosed  with  x-ray  mammography  or  a  combination  of  x-ray  mammography  and  other  imaging  modalities  (ultrasound, 
MRI)  and  biopsy.  The  patients  were  imaged  in  the  surgery  room  before  surgery.  Two-dimensional  optoacoustic  images 
were  acquired  in  several  locations  on  the  breast  with  cancer.  Irradiation  with  16  pulses  from  a  nanosecond  NdrYAG  laser 
operating  at  the  wavelength  of  1064  nm  was  performed  only  at  one  site  of  breast  skin  surface.  Optoacoustic  transducer  array 
was  placed  on  the  opposite  side  of  the  breast  approximately  beneath  the  point  of  irradiation.  Irradiation  point  was  placed 
approximately  above  the  area  suspicious  of  being  a  tumor.  The  exact  location  of  tumors  was  not  known  prior  to  opto¬ 
acoustic  imaging  procedure,  however,  approximate  location  could  be  determined  from  the  mammogram.  Tumors  were 
sometimes  palpable  and  biopsy  incision  was  visible  on  the  spared  segment  of  skin.  However,  some  tumors  were  not  palpable 
and  located  deep  within  the  breast.  Two  exemplary  optoacoustic  images  of  breast  carcinoma  are  presented  in  Figure  5. 
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Figure  5.  Optoacoustic  image  of  a  malignant  tumor  obtained  In  vivo  (  lett,  UTMB  patient  #  315221Q),  Optoacoustic  image  of  a  malignant 
tumor  obtained  in  vivo  (right,  UTMB  patient  #  026372P). 


Figure  6.  X-ray  mammography  image  of  a  cancerous  breast  (left).  There  were  2  suspicious  areas  on  mammogram,  only  one  of  them  was  a 
ductal  carcinoma.  Ultrasound  image  of  the  same  breast  showing  malignant  tumor.  UTMB  Patient  #315221 . 


The  main  conclusion  that  could  be  made  from  clinical  experiments  on  breast  cancer  patients  is  that  optoacoustic  tomography 
provides  enhanced  contrast  between  normal  tissues  and  cancerous  tumors.  This  contrast  varies,  however,  no  one  tumor  of  5 
tumors  examined  had  contrast  less  than  1,  i.e.  the  optoacoustic  amplitude  is  twice  higher  in  tumors  relative  to  normal 
background.  Optoacoustic  images  presented  in  Fig.  5  depict  details  of  tumor  structure  based  mainly  on  optical  properties  of 
cancerous  tissues.  Location  of  the  tumors  inside  the  breast  and  the  tumor  dimensions  also  were  accurately  determined  from 
the  optoacoustic  images  and  confirmed  through  comparison  with  x-ray  mammography,  ultrasound  and  pathological 
examination  (see  Figure  6).  Figure  6  presents  x-ray  mammography  and  ultrasound  images  of  the  breast  with  a  malignant 
tumor,  which  was  depicted  in  optoacoustic  image  in  Fig.  5  (left).  The  position  of  the  tumor  relative  to  the  skin  surface  was 
better  determined  from  ultrasound  image.  However,  detailed  structure  of  the  tumor  was  better  visualized  in  x-ray 
mammography  image.  On  the  other  hand,  x-ray  image  depicted  two  areas  suspicious  of  being  cancerous,  and  only  one  area 
contained  malignant  tumor  (as  determined  by  pathology  study).  The  tumor  core  brightly  displayed  in  the  optoacoustic  image 
was  not  visible  in  neither  x-ray  or  ultrasound  images,  but  was  confirmed  after  subsequent  pathology  examination.  The 
dimension  of  a  tumor  core  can  be  estimated  as  10x7  mm.  Tumor  was  detected  at  the  depth  of  11  mm  from  the  surface. 
Contrast  between  the  tumor  and  surrounding  normal  tissues  in  this  image  exceeds  1.2.  The  contrast  in  ultrasound  image  and 
x-ray  mammography  image  of  the  same  tumor  did  not  exceed  0.1.  Resolution  of  image  visualization  with  LOIS  (~  1-mm) 
was  comparable  with  that  of  x-ray  mammography  and  ultrasound.  Results  of  comparison  between  optoacoustic,  ultrasonic 
and  x-ray  radiography  images  for  the  patient  026372  (see  Fig.  5  ,  right)  and  for  the  other  3  patients  were  qualitatively  similar. 


3.  CONFOCAL  OPTOACOUSTIC  MICROSCOPY 


Confocal  Optoacoustic  transducer 

Confocal  optoacoustic  transducer  has  sharp  focusing  of  the  incident  optical  beam  and  narrow  focus  of  the  ultrasonic  detection 
system.  The  distribution  of  focused  light  and  the  caustic  of  the  ultrasonic  detection  form  the  volume  of  monitoring  for  the 
COAT.  The  schematic  diagram  of  COAT  is  presented  in  Fig.  1.  Pulsed  laser  radiation  was  delivered  to  the  COAT  via 
optical  fiber ,  focused  by  a  lens  onto  the  tissue  surface  through  opto~acoustic  lens.  Ultrasonic  waves  induced  in  tissue  by 
laser  pulses  propagating  backward  through  the  opto-acoustic  lens,  were  focused  with  the  opto-acoustic  lens  onto  the  ring- 
shaped  piezoelectric  detector.  Combination  of  flat  piezo- transducer  and  opto-acoustic  lens  is  aimed  to  obtain  focused 
receiving  ultrasonic  transducer  with  narrow  and  long  waist  of  the  focal  zone  right  under  the  irradiated  tissue  surface.  The 
shape  of  the  lens  curvature  allowed  transformation  of  a  spherical  acoustic  wave  from  tissue  into  a  planar  acoustic  wave 
without  phase  distortion.  A  charge  preamplifier  was  used  to  transform  acoustic  transducer  signals  into  amplified  electronic 
signals  for  further  data  acquisition  and  processing. 

Laser  pulses  from  a  Q-switched  Nd:YAG  laser  operating  at  the  fundamental  (A.=1.064  lam),  the  second  (1=0.532  |Lim)  and  the 
third  (1=0.355  fom)  harmonics  were  used  in  studies  of  light  absorption  distribution  in  phantoms  and  tissues.  The  temporal 
shape  of  laser  pulses  was  close  to  Gaussian  with  full  width  at  the  max/e  level  2TL=12ns  for  1=532  nm  and  2TL=10ns  for  1=355 
nm.  The  energy  of  each  laser  pulse  was  measured  with  a  calibrated  joulemeter  (ED-200,  GenTec,  Canada).  Sensitivity  of  the 
opto-acoustic  transducer,  permitted  detection  of  ultrasonic  waves  with  signal- to-noise  ratio,  SNR=1  in  tissues  with  absorption 
of  20  cm  ^  while  using  incident  laser  fluence  of  Fq<1  mJ/cm^,  suitable  for  safe  application  in  medical  imaging.^^ 
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Figure  7.  Principle  schematics  of  Confocal  Opto- Acoustic  Transducer. 

3,2.  Test  experinments  and  tissue  phantoms 

Opto-acoustic  imaging  studies  utilizing  confocal  optoacoustic  transducer  were  carried  out  with  optically  clear  solutions  of 
potassium  chromate  and  optically  turbid  (scattering)  aqueous  solutions  of  polystyrene  microspheres  colored  with  potassium 
chromate.  A  solution  of  3.5  g  of  potassium  chromate  in  10  cm^  of  water  yields  an  absorption  coefficient  of  about  10"*  cm'^  at 
the  wavelength  of  355  nm.  Optical  absorption  coefficient  of  aqueous  solutions  at  the  wavelength  of  355  nm  was  varied  in  the 
range  from  70  to  1,000  cm’^  in  order  to  study  response  of  the  opto-acoustic  transducer  to  ultrasonic  signals  with  various 
temporal  shape  and  duration.  Acoustic  properties  of  solutions  with  different  concentrations  of  potassium  chromate  were 
similar. 


3.2.1.  Spectral  sensitivity  of  COAT. 


The  spectral  response  of  ultrasonic  sensitivity  for  COAT  is  presented  in  Figure  8.  This  curve  was  obtained  irradiating  a  thin 
layer  of  concentrated  aqueous  solution  of  potassium  chromate  cm'^)  placed  onto  the  surface  of  the  Opto- Acoustic 

lens.  The  absorbing  solution  was  irradiated  with  a  wide  laser  beam.  The  spectral  response  of  COAT  was  calculated  as  a  ratio 
of  spectra  of  the  detected  opto-acoustic  signal  and  the  Fourier  spectrum  of  laser  pulse  measured  with  fast  photodiode.  The 
spectral  transfer  function  of  laser  generation  of  ultrasound  for  media  with  very  strong  optical  absorption  is  constant. 
Therefore,  laser- generated  spectrum  of  opto-acoustic  transients  should  coincide  with  the  spectrum  of  the  laser  pulse. 
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Figures.  Absolute  sensitivity  of  COAT 


3.1.2.  Reconstruction  of  light  absorption. 

A  homogeneous  absorbing  media  were  employed  to  demonstrate  the  procedure  of  reconstruction  of  absorbed  optical  energy 
distribution.  Aqueous  solution  of  potassium  chromate  was  irradiated  by  the  third  harmonic  of  NdrYAG  laser.  The  de- 
convolution  of  detected  signal  and  reference  OA-signal,  obtained  with  high  absorbing  medium  (light  absorption  coefficient 
about  10"^  cm'^),  was  applied.  Figure  9  shows  an  example  signal  detected  in  potassium  chromate  solution  with  absorption 
coefficient  of  (6.5±0.1)*10^  cm'^  (line  with  dots).  The  profile  of  the  optoacoustic  transient  depicted  in  Fig.  9  is  determined  by 
the  optical  absorption  coefficient,  the  laser  pulse  duration,  diffraction  of  ultrasonic  wave  and  transducer  sensitivity  within  the 
detected  frequency  band.  The  result  of  de-convolution  is  presented  on  Fig.  9  as  a  plain  line.  In  clear  absorbing  aquesous 
solutions,  the  exponential  slope  of  the  opto-acoustic  profile  is  defined  solely  by  the  optical  absorption  coefficient. 
Exponential  fitting  the  slope  yielded  the  absorption  coefficient,  which  was  equal  to  (648±8)  cm  ^ 


Figure  9.  Opto-acoustic  signal  (line  with  dots)  and  reconstructed  distribution  of  absorbed  optical  energy. 


3.2.3.  Image  of  cylindrical  objects. 


Polystyrene  latex  beads  with  average  diameter  of  0.76  jam  in  ~1%  concentration  were  employed  to  make  optically  scattering 
solutions  with  effective  scattering  coefficient  of  10  cm\  The  first  model  object  was  aluminum  wire  with  diameter  of 
approximately  lOO-pm  placed  in  optically  scattering  solution  at  various  depths  from  100  pm  to  1.5  mm.  The  second  model 
object  was  a  human  hair  embedded  in  aqueous  solution  at  the  depth  of  600-pm.  Thickness  of  the  hair  was  measured  with 
micrometer  and  equal  80  pm.  Thin  wire  and  hair  were  imaged  with  COAT  in  order  to  determine  in-depth  and  lateral 
resolution  of  the  confocal  opto-acoustic  front  surface  transducer. 

The  A1  wire  was  placed  in  polystyrene  beads  water  suspension  at  various  depths  parallel  to  the  transducer  surface.  Scans 
across  the  wire  axis  were  performed.  Opto-acoustic  signal  amplitude  as  a  function  of  lateral  coordinate  is  given  in  Figure  10. 
The  full  width  at  1/e  of  the  amplitude  equals  to  210  pm  for  290  pm  depth,  as  measured  from  the  curve  in  Fig.  10.  The 
dependence  widened  with  depth  of  detection.  Therefore  lateral  resolution  decreases  with  depth,  as  the  opto-acoustic  cone  of 
detection  widens. 

The  axial  in-depth  resolution  within  the  transducer  sensitive  volume  is  equal  to  the  thickness  of  tissue  where  laser-induced 
pressure  is  confined  during  laser  energy  deposition.  This  thickness  is  defined  by  the  product  of  the  laser  pulse  duration  and 
the  speed  of  sound  in  tissue  and  equal  18-pm. 
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Figure  10.  Signal  amplitude  dependences  vs  lateral  displacement  for  various  depth. 

To  demonstrate  spatial  resolution  of  the  transducer  opto-acoustic  imaging  of  human  hair  was  carried  out.  Hair  with  diameter 
of  80  pm  was  placed  in  water  and  irradiated  with  green  laser  pulses.  Scanning  across  the  hair  axis  was  performed  with  lateral 
step  of  25.4  pm.  Obtained  signals  were  processed  using  absorption  reconstruction  procedure  and  resulting  image  was  created. 
The  image  is  presented  in  Figure  1 1 .  It  shows  the  hair  as  a  bright  strip  with  the  thickness  corresponding  to  the  thickness  of 
the  hair,  as  in-depth  resolution  is  high.  The  width  of  the  image  is  several  times  greater,  then  the  thickness  due  to  relatively 
low  angular  aperture  of  the  transducer.  One  can  conclude,  that  in  confocal  laser  optoacoustic  microscopy  the  size  of  image  in 
lateral  direction  will  be  greater  then  in-depth  size  according  to  angular  aperture  of  optoacoustic  lens. 
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Figure  1 1 .  Opto> acoustic  image  of  human  hair. 

3.3.  Nonin vasive  detection  of  oral  cancer  in  vivo 


A  well-established  animal  model  of  DMBA-induced  squamous  cell  carcinoma  in  oral  cavity  was  used  in  this  study  [35]. 
Syrian  golden  hamsters,  weighing  -200  grams  were  employed.  The  motivation  of  using  hamster  model  was  the  development 
of  hamster  buccal  pouch  carcinoma  comparable  to  human  oral  cancer,  as  shown  in  [36,37].  DMBA  solution  (0.5%  of  9, 10- 
Dimethyl- 1,2-Benzanthracene  dissolved  in  mineral  oil)  was  applied  twice  a  week  onto  the  buccal  mucosa  of  the  left  cheek 
pouch  using  a  small  brush  until  the  first  lesion  was  visible  with  a  naked  eye.  Photographs  of  the  final  lesions  were  taken  and 
the  animal  was  euthanized.  Hamster  pouches  were  cut,  stretched  and  fixed  with  pins  on  a  flat  cork,  and  then  embedded  in 
paraffin.  Tissue  sections  were  stained  with  H&E  (hematoxylin-eosin).  Histological  sections  were  used  as  a  gold  standard  in 
order  to  compare  with  opto-acoustic  images  of  different  types  of  lesions  and  stages  of  squamous  cell  carcimona  (from 
displasia  to  carcinoma  in  situ  to  invasive  malignant  lesion). 

Opto-acoustic  imaging  scans  were  performed  every  week  during  the  lesion  induction  period.  The  procedure  consisted  of 
translating  the  transducer  along  the  wet  surface  of  pouch  mucosa.  A  slight  pressure  was  applied  for  better  acoustic  contact 
that  resulted  in  flattening  of  irradiated  tissue  surface.  Confocal  opto-acoustic  transducer  was  applied  to  the  pouch  surface 
with  a  thin  water  layer  to  provide  acoustic  coupling  between  the  transducer  and  the  tissue.  During  the  experimental  period, 
the  animals  had  free  access  to  food  and  water,  except  during  the  actual  opto-acoustic  measurements,  when  no  food  was 
allowed  in  order  to  prevent  potential  contamination  of  the  opto-acoustic  images. 

The  main  goal  of  our  study  was  to  test  feasibility  of  the  proposed  technique  in  detection  and  staging  of  cancer  in  oral  mucosa. 
Previous  studies  performed  with  well-defined  models  of  animal  oral  carcinoma  consistently  reported  the  following  sequence 
of  histologic  changes  that  occur  in  the  affected  pouch  mucosa:  (1)  hyperplasia,  (2)  displasia,  (3)  carcinoma  in  situ,  (4) 
invasive  carcinoma  and  (5)  extensive  tumors.  Scanning  of  the  opto-acoustic  transducer  along  the  cheek  surface  was 
performed  to  acquire  data.  Images  presented  below  show  a  thin  (10-30  pm)  water  layer  as  the  top  layer,  so  that  the  first 
mucosal  surface  is  located  at  about  10-30-pm  depth.  The  recorded  signals  were  processed  using  reconstruction  procedure 
(see  section  2.4)  and  plotted  on  the  two  dimensional  plane  as  25-pm  wide  columns.  Thereby  reconstruction  of  two- 
dimensional  images  was  made  in  hamsters  with  various  stages  of  cancer  development.  Opto-acoustic  images  were  compared 
with  histological  microscopic  sections  stained  with  hematoxylin-eosin. 

Optoacoustic  imaging  of  normal  mucosa  in  oral  cavity 

Histological  H&E  sections  of  normal  hamster  cheek  pouch  were  described  in  the  literature  previously. There  are  several 
distinct  layers  in  a  normal  cheek  pouch.  Figufe  12  depicts  a  microscopic  histology  section  (left)  and  corresponding  opto- 
acoustic  image  of  a  normal  hamster  cheek  pouch  (right).  Both,  the  histology  and  the  opto-acoustic  tomography  display 
distinct  layers  that  correspond  to  keratinized  stratified  squamous  epithelium  (E),  a  very  thin  basal  membrane  separating 
epithelium  from  mucosa  of  lamina  propria  (L),  muscle  fiber  mucosa  (M),  submucosa  connective  tissue  (S). 
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Fig.  12.  Histological  section  and  opto-acoustic  image  ot  hamster  normal  cheek  pouch.  Normal  buccal  mucosa  pouch  is  characterized 
layered  structure  with  distinct  layers:  keratin  on  top  of  epithelium  (E),  lamina  propria  (L)  muscle  fiber  mucosa  (M),  submucosa  connective 
tissue  (S). 


3.4.2  Displasia 

A  process  of  cancer  development  leads  to  gradual  destruction  in  layered  structure  of  mucosa.  The  earlies  stage  of  cancer  is 
characterized  by  hyperplasia  and  irregularly  increased  thickness  of  epithelium,  however  without  compromise  to  the  layered 
structure.  Hyperplasia  is  an  early  precancerous  stage,  which  can  not  be  defineteiy  outlined. 
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Figure  13.  Earlies  stage  of  cancer  development.  An  opto-acoustic  image  measured  in  vivo  from  a  normal  hamster  cheek  pouch  (right 
image)  and  corresponding  microscopic  histology  section  (left  image)  employing  hematoxylin-eosin  staining. 

The  second  stage,  displasia,  is  characterizaed  with  more  prominent  thickening  of  epithelium  layer  with  rare  and  irregular 
breaks  in  the  layered  structure,  made  of  local  accumulation  of  epithelial  cells  with  bigger  nuclei  and  higher  concentration  of 
nuclei.  Histological  criteria  used  for  determining  dysplasia  is  irregular  epithelial  stratification  These  abnormal  epithelial  cells 


start  proliferation  into  deeper  mucosal  layers.  Histological  section  of  dysplasia  is  depicted  in  Figure  8  (left).  Opto-acoustic 
image  presented  in  Fig.  13  (left)  shows  a  focal  point  of  displasia.  The  epithelium  layer  looks  almost  normal  around  the 
microscopic  region  with  displasia. 


3.4.3.  Carcinoma  in  situ 

Figure  14  presents  histological  section  and  opto-acoustic  image  of  hamster  pouch  treated  with  DMBA  for  6weeks  (left 
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Figure  14.  A  microscopic  histology  section  showing  carcinoma  in  siiu  ol  a  hamster  cheek  pouch  treated  with  DMBA  carcinogenic  agent 
(left  image,  H&E  staining)  and  corresponding  opto-acoustic  image  (right  image). 

Microscopic  squamous  cell  carcinoma  in  situ  were  noted  in  hystology  section  by  presence  of  acanthosis,  lack  of  maturation, 
cellular  atypia,  mitoses.  Opto-acoustic  image  (right)  shows  partial  loss  of  layered  tissue  structure  with  local  invasions  of 
epithelial  ceils  into  deeper  mucosa,  irregularly  shaped  areas  in  mucosa  of  lamina  propria. 
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Figure  15.  H&E  stain  of  DMBA  model  treated  during  12  weeks  (left)  and  corresponding  opto-acoustic  image  (right).  Histological  section 
of  the  hamster  pouch  shows  invasive  carcinoma.  The  tumor  nests  have  an  irregular  outline  with  infiltrative  borders. 


Further  progression  of  carcinoma  development  was  observed  in  all  animals  by  presence  of  invasive  carcinoma.  Hamster 
pouch  treated  with  DMBA  for  12  weeks  showed  structural  changes  characterized  as  multiple  carcinomas  in  situ  and  invasive 
carcinomas.  Scanning  across  lesion  area  was  made  and  opto-acoustical  image  was  reconstructed.  A  typical  opto-acoustic 
image  showing  invasive  carcinoma  and  correspondent  histology  H&E  section  are  displayed  in  Figure  15). 

CONCLUSIONS 

The  first  clinical  prototype  of  the  laser  optoacoustic  system  for  two-dimensional  imaging  of  the  breast  in  vivo  was  developed 
and  employed  for  detection  of  breast  cancer  in  5  patients.  The  system  employed  an  array  of  32  elements  of  ultra  wide-band 
PVDF  transducers  (20kHz-2MHz)  shaped  as  an  arc  of  120  degrees.  The  system  initially  was  tested  in  various  phantoms  and 
demonstrated  specifications  that  make  it  suitable  for  high-contrast  imaging  human  tissues.  Sensitivity  of  LOIS  permits 
detection  of  2-mm  blood  vessels  at  the  depth  of  7.5  cm.  In  depth  resolution  equals  0.4  mm,  lateral  resolution  is  about  1  mm 
depending  on  position  of  tumor  relative  to  the  transducer  array.  Quantitative  information  on  distribution  of  absorbed  optical 
energy  can  be  obtained  from  optoacoustic  images.  This  clinical  prototype  LOIS  is  the  basis  for  the  next  generations  of  three- 
dimensional  imaging  systems  to  be  developed  at  UTMB  in  collaboration  with  LaserSonix  Technologies,  Inc.  (Galveston, 
TX). 

A  novel  confocal  opto-acoustic  transducer  for  imaging  layered  tissues  in  vivo  was  developed.  This  transducer  combines 
focused  optical  system  for  laser  irradiation  of  tissue  and  a  focused  piezoelectric  detection  system  operating  in  the  ultrasonic 
band  of  IMHz  tolOO  MHz. 

Axial  in-depth  resolution  of  the  transducer  is  15  jim  for  biological  tissue.  Lateral  resolution  is  not  worse  than  60  pm  at  the 
depth  of  up  to  600  pm.  The  major  advantage  of  COAT  compared  with  OAFST  is  improved  lateral  resolution  over  the  range 
of  depths  in  tissue  from  0  to  1.5  mm.  A  deconvolution  procedure  for  reconstruction  of  absorbed  energy  distribution  was 
developed  and  tested  yielding  satisfactory  accuracy  of  quantitative  analysis. 

Confocal  opto-acoustic  imaging  was  sufficiently  sensitive  to  depict  structural  changes  in  mucosal  tissues  much  earlier  than 
any  changes  could  be  noticed  on  the  tissue  surface  with  naked  eye.  Furthermore,  the  opto-acoustic  images  yielded 
information  of  in-depth  tissue  structure  that  closely  resembled  structures  provided  by  optical  microscopy  of  histology 
sections.  Microscopic  displasia  can  be  differentiated  from  the  normal  mucosa  and  carcinoma  in  situ  could  be  differentiated 
from  displasia.  The  differentiation  of  the  invasive  carcinoma  from  normal  and  precancerous  lesions  was  very  prominent. 
Changes  associated  with  gradual  cancer  development  could  be  characterized  as  a  transformation  of  tissue  layered  structure 
into  a  heterogeneous  structure.  The  nature  of  the  opto-acoustic  contrast  in  mucosa  and  carcinoma  needs  to  be  established  in 
order  to  enhance  further  the  imaging  sensitivity.  We  speculate  that  differentiation  of  cancer  is  based  on  tumor  angiogenesis. 

The  results  of  the  presented  studies  confirmed  feasibility  of  COAT  technique  in  differentiation  of  normal  mucosa  from 
cancerous,  and  in  characterization  of  early  stages  of  squamous  cell  carcinoma  in  animal  model.  The  results  from  the  in  vivo 
study  suggest  that  there  is  a  potential  for  confocal  opto-acoustic  tomography  for  imaging,  monitoring  and  guiding  biopsy. 
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